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ABSTRACT 

This  thesis  developed  a  data  acquisition  system  to  be 
used  in  conjunction  with  the  3.5'  x  5.0'  low  speed  wind 
tunnel  at  the  Naval  Postgraduate  School.  Interactive 
graphic  programs  were  developed  to  aid  in  data  acquisiton 
and  analysis.  In  addition,  the  internal  balance  that  was 
designed  by  Major  Scott  Mair  and  Major  Chris  Sargent  was 
redesigned  to  correct  some  problems  encountered  with  the 
drag  component.  The  balance  was  also  instrumented  to  record 
the  pitch  and  yaw  moment  components.  A  calibration  rig  was 
designed  and  constructed  in  order  to  evaluate  the 
interactions  of  the  different  components.  The  equipment 
used  and  programs  developed  for  data  acquisition  and 
analysis  were  adequate.  However,  balance  calibration 
revealed  problems  with  the  calibration  rig  and  location  of 
the  roll  component  strain  gage.  Both  of  these  problems  will 
have  to  be  corrected  before  accurate  readings  can  be 
expected  from  this  balance  design. 
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I.   INTRODUCTION 

A.   BACKGROUND 

This  project  used  the  3.5'  x  5.0'  low  speed  wind  tunnel 
at  the  Naval  Postgraduate  School  to  continue  the 
aerodynamic  study  of  the  effects  of  helicopter  noses  and 
tails  on  drag  conducted  by  Major  Mair  [Ref.  1]  and  Major 
Sargent  [Ref.  2].  Their  work  included  the  construction  of 
nine  various  helicopter  configurations,  a  sting  mounted 
support  system  and  a  four-degree-of-freedom  balance  that 
was  internal  to  the  model.  They  also  studied  the  airflow 
around  the  various  configurations  with  cotton  tufting. 

To  reduce  and  analyze  the  data,  Majors  Mair  and  Sargent 
developed  several  independent  computer  programs.  These 
programs  calibrated  the  balance,  corrected  the  raw  data 
from  the  wind  tunnel  runs,  computed  the  equivalent  flat 
plate  area  and  produced  plots  of  various  parameters  for 
comparison.  However,  due  to  problems  encountered  with  the 
drag  component  of  the  balance,  their  results  were 
inconclusive. 

B .   GOALS 

The  primary  goal  of  this  project  was  to  provide  the 
students  enrolled  in  the  helicopter  design  class  at  the 
Naval  Postgraduate  School  with  a  laboratory  type  set-up 
that  would  develop  realistic   Equivalent   Flat   Plate  Area 


information  for  various  helicopter  configurations.  This 
included  creating  interactive  computer  programs  that  the 
students  could  use  with  the  personal  computer  located  at 
the  wind  tunnel. 

The  secondary  goals  of  this  project  were  to  upgrade  the 
two-axis  internal  wind  tunnel  balance  to  a  six-axis 
internal  balance  and  to  develop  a  calibration  rig  to 
evaluate  the  interactions  of  the  balance  components. 

Three  landing  gear  configurations  were  also  designed 
and  constructed  for  future  analysis  of  the  drag  that  they 
add  to  the  helicopter. 
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II.  APPROACH  TO  THE  PROBLEM 

A.   LANDING  GEAR  DESIGN 

To  provide  a  realistic  representation  of  landing  gear 
used  with  modern-day  helicopters,  one  type  of  landing  gear 
was  selected  for  each  nose  shape  (Figures  2.1-2.3). 
Planviews  for  the  three  types  of  landing  gear  were  prepared 
and  are  included  in  Appendix  A. 

For  the  attack  nose,  a  skid  type  of  landing  gear  was 
constructed  of  aluminum  tubing  secured  to  an  aluminum 
plate.  This  type  of  landing  gear  is  considered  a  fixed  gear 
but  was  chosen  because  of  its  wide  use  for  numerous 
helicopters.  For  both  the  smooth  nose  and  blunt  nose  a 
simulated  retractable  gear  was  constructed  of  wood 
stubwings  and  model  airplane  tires.  Threaded  inserts  were 
mounted  in  the  noses  and  stubwings  to  allow  easy  removal  of 
the  wheel  and  strut  assemblies.  This  allowed  the  models  to 
be  tested  in  both  the  clean  and  dirty  configuration. 

B.   INTERNAL  BALANCE  MODIFICATION 

A  modified  Mair/Sargent  balance,  Figures  2.4-2.7,  was 
used  for  this  project. 

To  improve  the  output  recorded  from  the  axial 
component,  the  cuts  (see  Figure  2.4)  alongside  that  cavity 
were  increased  by  1/16  of  an  inch.  In  addition,  the  cavity 
itself   was   squared  off  thus  reducing  the  curvature  of  the 

11 


Figure  2.1   Attack  Nose  with  gear 


f«» 


Figure  2.2    Smooth  Nose  with  gear 
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Figure  2.3    Blunt  Nose  with  gear 
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Figure  2.4    Internal  Balance  with  Modifications 


Figure  2.5    Re-Wired  Internal  Balance 
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surfaces  in  the  cavity.  These  modifications  increased  the 
flat  area  upon  which  the  strain  gages  could  be  placed. 

An  area  (see  Figure  2.6),  similar  to  that  for  the 
pitching  moment  component,  was  cut  to  allow  recording  of 
the  yawing  moment  component. 

For  compatability  and  increased  accuracy,  the  aluminum 
strain  gages  were  removed  and  replaced  by  EA-09-062AQ-350 
stainless  steel  ones.  •  These  gages  were  smaller  and  thus 
allowed  for  better  placement  within  the  cavities.  To 
provide  a  longer  life,  the  gages  were  bonded  to  the  balance 
with  an  M-Bond  AE-15  adhesive  system.  They  were  cured  at  a 
temperature  of  150  degrees  Fahrenheit  for  two  hours. 

It  was  desired  to  record  the  six  component  forces  on 
the  helicopter;  lift,  drag,  yaw,  pitching  moment,  yawing 
moment  and  rolling  moment.  However,  since  the  internal 
balance  was  designed  to  record  only  four  components,  the 
sting  support  was  instrumented  to  record  the  yaw  force  and 
rolling  moment  (Figure  2.8). 

The  gages  for  the  yaw  component  were  placed  on  the  side 
of  the  sting  support  to  undergo  tension  and  compression 
when  subjected  to  a  yawing  force.  The  gages  for  the  rolling 
moment  component  were  placed  on  top  and  bottom  of  the  sting 
support  at  a  45  degree  angle  to  the  sting  axis.  Thus,  they 
experienced  tension  and  compression  when  the  model  was 
subjected  to  a  rolling  moment. 
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Figure  2.6    Lift,  Pitch  Moment  and  Yaw  Moment  Gages 


Figure  2 . 7    Drag  Component  Gage 


Yaw  and  Roll  Moment  Gages 
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C.   BALANCE  CALIBRATION 

In  order  to  determine  any  interaction  between  the  six 
recording  components,  a  calibration  rig  was  designed  and 
mounted  to  the  main  fuselage  (Figure  2.9).  The  pans  and 
pullies  were  arranged  to  simulate  the  twelve  forces  and 
moments  that  the  helicopter  models  would  experience  in  the 
wind  tunnel.  Each  pan,  or  component  was  loaded  from  zero  to 
twenty  pounds  in  one  pound  increments.  With  each  loading, 
readings  from  the  other  five  components  were  recorded.  This 
method  was  conducted  until  all  six  components  were  loaded 
in  both  the  positive  and  negative  direction  and  produced 
thirty  interaction  matrices. 

A  balance  calibration  program,  Figure  A. 16,  was  written 
to  determine  the  relations  required  to  convert  the  raw  data 
counts  to  actual  forces  and  moments.  For  each  component 
loading,  the  prime  gage  constants  were  determined  using  the 
following  least  squares  curve  fit  [Ref.  3]: 

[Sum(Xi) A2]*K1  +  [Sum(XiA3) ] *K2  =  Sum(Xi*Yi) 

[Sum(Xi) A3]*K1  +  [Sum(XiA4) ] *K2  =  Sum ( (XiA2 ) *Yi)  . 

where  Xi  is  the  raw  data  count 
and  Yi  is  the  applied  load 

Once  these  constants  were  computed  for  all  twelve  loadings, 
the  raw  data  counts  were  converted  to  forces  and  moments. 
Then  for  each  interaction  matrix,  the  same  least  squares 
equations  were  used  to  determine  the  interaction 
coefficients.   These   prime   gage  constants  and  interaction 


16 


Figure  2.9    Balance  Calibration  Rig 
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coefficients  were  stored  for  use  in  the  reduction  of  the 
raw  data  recorded  during  the  wind  tunnel  tests. 

D.   DATA  ACQUISITION 

To  improve  data  reduction,  it  was  desired  to  have  a 
data  acquisition  system  built  into  the  wind  tunnel 
instrumentation.  The  strain  gages  on  the  internal  balance 
and  support  sting  were  connected  to  a  Pacific  Instruments 
model  8255  transducer  amplifier  (Figure  2.10).  The  outputs 
from  this  unit  were  then  routed  to  data  acquisition  cards 
mounted  in  an  IBM  PC  AT.  A  data  acquisition  program,  Figure 
A. 10,  was  incorporated  into  the  interactive  program  set 
that  allowed  for  quick  and  easy  recording  of  test  data. 

Having  the  transducer  amplifier  connected  to  the  strain 
gages  allowed  for  easy  zeroing  and  balancing  of  the  gages. 
The  data  acquisition  set-up  provided  real  time  readout  of 
the  forces  that  the  helicopter  was  experiencing  in  the  wind 
tunnel . 
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Figure  2.10   Data  Acquisition  Equipment' 
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III.  SOLUTION  TO  THE  PROBLEM 

A.   DATA  COLLECTION  AND  REDUCTION 

The  data  collection  program,  Figure  A. 10,  was  designed 
to  record  the  data  in  a  standardized  collection  method. 
Interactive  steps  instruct  the  user  to  follow  the 
procedures  listed  below. 


1) .   Zero  all  amplifiers  without  the  model  in  place. 
2 ) .   Mount  the  model  to  the  sting  support  and  record 

a  zero  force  reading. 
3).   Place  the  calibration  switches  to  the  +  position 

and  record  a  calibration  reading. 
4) .   Replace  the  calibration  switchres  to  the  center 

position  and  start  the  wind  tunnel. 
5) .   Record  first  and  last  data  point  at  zero  angle  of 

attack  with  all  other  readings  between  plus  ten 

degrees  and  minus  eight  degrees  angle  of  attack. 


By  using  a  delay  loop,  changes  could  be  made  to 'the 
angle  of  attack  of  the  model  before  the  next  set  of  data 
points  were  recorded.  For  each  tunnel  speed,  the  data 
recorded  for  the  various  angles  of  attack  were  stored  on 
separate  files. 

Once  the  tunnel  runs  were  completed  and  the  data  was 
recorded  on  file,  another  program,  Figure  A. 11,  was 
designed  to  convert  the  raw  readings  to  readings  of  forces 
and  monemts.  The  sixty  interaction  coefficients  developed 
from  the  balance  calibration  were  incorporated  into  this 
data  reduction.   The  following  equations,  obtained  from  Mr. 
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David  Backs  at  the  NASA  AMES  Reseasrch  Center,  were  used  to 
correct  for  the  interaction  of  the  balance  components. 

Fl'  =  Fl  -  (delFl/delF2)*F2  -  (delFl/delF3) *F3  -  

-  (delFl/delF2A2)*F2A2  -  - (delFl/delF6A2) *F6A2 

Fl • ' =  Fl  -  (delFl/delF2 ) *F2 • - 

-  (delFl/delF2A2) *F2,A2  - (delFl/delF6A2) *F6 ' A2 

These  equation  were  written  for  each  component  and 
placed  in  an  iterative  loop  that  checked  the  difference 
between  the  two  prime  values.  After  the  forces  were 
corrected  for  component  interaction,  a  weight  tare  equation 
was  used  to  correct  for  the  weight  of  the  model 

B.   ANALYSIS  OF  DATA 

The  main  file  of  the  interactive  programs,  Figure  A. 9, 
controlled  the  data  recording,  data  reduction  and  data 
analysis  with  a  menu  format.  For  data  analysis,  the  stored 
converted  data  readings  were  used  to  calculate  the  lift  and 
drag  coefficients  and  equivalent  flat  plate  area  for  each 
angle  of  attack.  Files  were  created  for  coefficient  of  lift 
versus  angle  of  attack,  coefficient  of  drag  versus 
coefficient  of  lift,  coefficient  of  drag  versus  coefficient 
of  lift  squared  and  equivalent  flat  plate  area  versus  angle 
of  attack.  There  was  one  file  of  each  created  for  each 
tunnel  speed. 

A  plotting  routine,  Figure  A. 15,  was  included  in  the 
interactive  programs  to  allow  quick  analysis  of  the 
recorded  data.   To  aid  in  the  analysis,   up  to  three   plots 
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could  be  shown  on  one  graph.  A  delay  loop  was  incorporated 
into  the  plotting  routine  to  allow  for  the  option  of 
obtaining  a  hard  copy  of  the  graph  by  using  the  print 
screen  command. 
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IV.    RESULTS 

A.   BALANCE  CALIBRATION 

The  balance  was  loaded  in  the  twelve  component 
directions  and  the  balance  calibration  program  was  used  to 
produce  the  tables  B.2  through  B.13.  During  the  balance 
calibration,  a  drift  in  each  of  the  components  was  noted 
under  steady  state  conditions.  To  correct  for  this  problem 
it  was  assumed  that  each  component  had  a  constant  drift 
rate.  The  difference  between  the  first  and  last  zero 
reading  was  divided  by  the  number  of  data  points  taken. 
This  correction  factor  was  then  applied  to  each  data  point. 
This  method  produced  good  correlation  between  different 
data  points  taken  for  the  same  load. 

B.   DATA  COLLECTION 

Test  runs  were  conducted  to  evaluate  the  data 
acquisition  program.  The  vibrations  of  the  model  in  the 
wind  tunnel  resulted  in  erratic  fluctuations  in  the  sense 
indicators  of  the  amplifiers.  To  correct  for  this  problem, 
the  acquisition  program  was  modified  to  collect  one  hundred 
samples  at  a  rate  of  five  hundred  samples  per  second.  The 
average  of  these  one  hundred  samples  was  taken  as  one  data 
point.  This  method  produced  constant  readings  for  data 
points  of  similar  conditions. 
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C.   DATA  REDUCTION  AND  ANALSYIS 

The  data  reduction  program  was  used  to  convert  the  raw 
test  data.  The  interaction  equations  diverged  instead  of 
converging  to  a  single  value.  This  pointed  out  a  problem 
with  the  interaction  coefficients.  Examination  of  the 
calibration  tables  revealed  extremely  large  interactions 
between  the  loading  of  the  lift  component  and  reaction  in 
the  roll  moment  component.  There  was  also  noted  a  large 
interaction  between  the  loading  of  the  pitch  moment 
component  and  reaction  in  the  lift  component. 

The  reduction  program  was  modified  to  correct  for  the 
interactions  between  the  lift  and  drag  components  only.  The 
reduction  program  was  again  executed  using  the  recorded 
test  data.  This  time  the  interactions  converged.  This 
proved  the  validity  of  the  interaction  equations  and 
confirmed  the  problems  with  the  balance  calibration. 

A  test  file  was  created  to  display  the  plotting 
capabilities  of  the  interactive  programs.  Figure  4.1  shows 
the  results. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.  CONCLUSIONS 

The  data  acquisition  system  and  programs  developed  were 
adequate  for  recording  the  forces  and  moments  experienced 
in  the  wind  tunnel.  The  extremely  large  interactions 
between  the  lift  and  roll  moment  components  and  between  the 
pitch  moment  and  lift  components  (see  Tables  B.2  -  B.13) 
precluded  any  analysis  of  wind  tunnel  data.  However,  the 
data  acquisition  system  and  balance  calibration  program 
warrant  further  research  and  development  of  the  balance 
system. 

The  interactive  programs  greatly  reduced  the  workload 
required  in  the  data  acquisition  and  analysis  phases.  They 
provided  a  quick  and  easy  means  for  the  analysis  of  the 
recorded  data. 

B.  RECOMMENDATIONS 

The  following  are  given  as  recommendations  to  improve 
the  calibration  of  the  internal  balance  system. 

1.   Balance  Modification 

The  location  of  the  roll  moment  strain  gage  was 
decided  to  be  the  cause  of  the  large  interaction  between 
the  lift  and  roll  moment  component.  This  gage  was  located 
on  the  back  end  of   the   sting   support   system.   The   roll 
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moment  component  was  calibrated  by  applying  a  torque  to  the 
center  section  of  the  model  that  was  mounted  to  the 
internal  balance.  When  the  lift  component  was  loaded,  a 
large  bending  moment  was  felt  by  the  roll  moment  gage  and 
thus  producing  the  large  interactions. 

By  placing  the  roll  moment  strain  gage  at  a 
forty-five  degree  angle  on  the  same  cut-out  section  as  the 
pitch  moment  gage,  the  torque  applied  to  the  center  section 
can  be  used  for  calibration.  Also,  the  loading  of  the  lift 
component  will  not  greatly  affect  the  reaction  of  the  roll 
moment  component. 

2 .   Calibration  Rig  Modification 

By  using  the  center  section  of  the  model  to  mount 
the  calibration  pans  for  the  pitch  moment  component,  the 
strains  were  incorrectly  transmitted  to  the  lift  component. 
A  separate  calibration  set-up  will  have  to  be  designed  that 
isolates  the  force  appplied  to  the  pitch  moment  component 
from  the  rest  of  the  balance. 
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APPENDIX  A 
SKETCHES  AND  PROGRAMS 


Figure  A.l   Nose  Gear  for  Smooth  Nose 
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Figure  A. 2    Main  Gear  for  Smooth  Nose 
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Figure  A. 3    Stubwing  for  Smooth  Nose  "(Top  View) 
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Figure  A. 4    Skid  Gear  and  Wing  for  Attack  Nose 
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Figure  A. 5    Skid  Gear  for  Attack  Nose   (Top  View) 


32 


Figure  A. 6   Nose  Gear  for  Blunt  Nose 
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Figure  A. 7    Main  Gear  for  Blunt  Nose 
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Figure  A. 8   Stubwing  for  Blunt  Nose  (Top  View) 
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10    REM MAIN.BAS  (MAIN  CONTROL  PROGRAM) — 

2  0    '  PATRICK  A.  WITT 
30    '  20  JULY  1985 

40    ' ■ ■ 

50    REM 
280   KEY  OFF 
29  0   SCREEN  1 

3  00   COLOR  1 

310   FOR  I  =  1  TO  2 

320   IF  I  =  1  THEN  FILL  =  11  ELSE  FILL  =  12 

33  0   IF  I  =  1  THEN  EDGE  =  9  ELSE  EDGE  =  10 

340   IF  I  =  1  THEN  Xl=106:  X2=0:  X3=56  ELSE  Xl=108 :  X2=2 : 

X3=58 
350   IF  I  =  1  THEN  Yl=48:  Y2=94:  Y3=140  ELSE  Yl=46:  Y2=92: 

Y3=138 
360   »  MH" 

370   DRAW  "C=EDGE;BM=X1;,=Y1;U3  0R4D12R12U12R4D3  0L4U12" 
380   DRAW  "L12D12L4BE1P=FILL; ,=EDGE;" 
390   ■  "E" 
400   Xl=  XI  +  26 
410   DRAW  "C=EDGE;BM=X1; ,=Y1;U30R20D4L16D9R8D4L8D9R16D4L2 

0BE1P=FILL; ,=EDGE;" 
420   '  "L" 
430   Xl=  XI  +  26 
440   DRAW  "C=EDGE;BM=X1; , =Y1 ;U30R4D26R16D4L20BE1P=FILL; , 

=EDGE ; " 
450   '  "O" 
460   Xl=  XI  +  26 
470   DRAW  "C=EDGE;BM=X1; ,=Y1;U30R20D30L20BR4BU4U22R12D22L 

12BG1P=FILL; , =EDGE ; " 
48  0   'NEXT  LINE 
490   '  "C" 
500   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R20D4L16D22R16D4L20BE1P= 

FILL; ,=EDGE;" 
510   ■  "O" 
520   X2=  X2  +  26 
530   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R20D30L20BR4BU4U22R12D22L 

12BG1P=FILL; ,=EDGE;" 
540   '  "N" 
550   X2=  X2  +  26 
560   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R4F12U12R4D30L4H12D12L4BE 

1P=FILL; ,=EDGE;" 
570   '  "F" 
580   X2=  X2  +  26 
590   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R20D4L16D6R8D4L8D16L4BE1P= 

FILL;  ,=EDGE;" 
600   '  "I" 
610   X2=  X2  +  26 

620   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R4D30L4BE1P=FILL; ,=EDGE;" 
630   ■  "G" 

Figure  A.  9   MAIN.BAS  -  Main  Controlling  Program 
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640   X2  =  X2  +  14 

650   DRAW  "C=EDGE;BM=X2; ,=Y2;U3  0R20D4L16D22R12U4L4U4R8" 

660   DRAW  "D12L2  0BE1P=FILL; , =EDGE ; " 

670   '  "U" 

680   X2=  X2  +  26 

69  0   DRAW  "C=EDGE;BM=X2; ,=Y2;U3  0R4D2  6R12U2  6R4D3  0L2  0BE1P= 

FILL;,=EDGE;" 
700   '  "R" 
710   X2=  X2  +  26 
720   DRAW  "C=EDGE;BM=X2; , =Y2 ;U30R20D15L12F15L4H15D15L4BU1 

9BR4U7" 
730   DRAW  "R12D7L12BL1P=FILL; ,=EDGE;" 
740   ■  "A" 
750   X2=  X2  +  30 
760   DRAW  "C=EDGE;BM=X2; , =Y2 ;U30R20D30L4U15L12D15L4BU19BR 

4U7" 
770   DRAW  "R12D7L12BL1P=FILL;,=EDGE;" 
780   •  "T" 
790   X2=  X2  +  26 
800   DRAW  "C=EDGE;BM=X2; ,=Y2;BU26U4R20D4L8D26L4U26L8BE1P= 

FILL;, -EDGE;" 
810   ■  "I" 
820   X2=  X2  +  26 

83  0   DRAW  "C=EDGE;BM=X2; ,=Y2;U3  0R4D3  0L4BE1P=FILL; ,=EDGE;" 

84  0   '  "0" 

850   X2=  X2  +  14 

860   DRAW  "C=EDGE;BM=X2; ,=Y2;U30R20D30L20BR4BU4U22R12D22L 

12BG1P=FILL; ,=EDGE;" 
870   •  "N" 
880   X2=  X2  +  26 
890   DRAW  "C=EDGE;BM=X2; ,=Y2;U3  0R4F12U12R4D3  0L4H12D12L4BE 

1P=FILL; ,=EDGE;" 
9  00   'NEW  LINE 
910   '  "A" 
920   DRAW  "C=EDGE;BM=X3; ,=Y3;U30R20D30L4U15L12D15L4BU19BR 

4U7" 
93  0   DRAW  "R12D7L12BL1P=FILL; , =EDGE ; " 
940   '  "N" 
950   X3=  X3  +  26 
960   DRAW  "C=EDGE;BM=X3; /=Y3;U30R4F12U12R4D30L4H12D12L4BE 

1P=FILL; ,=EDGE;" 
970   »  "A" 
980   X3=  X3  +  26 
990   DRAW  "C=EDGE;BM=X3; ,=Y3;U30R20D30L4U15L12D15L4BU19BR 

4U7" 
1000  DRAW  "R12D7L12BL1P=FILL; ,=EDGE;M 
1010  '  "L" 
1020  X3=  X3  +  26 
1030  DRAW  "C=EDGE;BM=X3; , =Y3 ;U30R4D26R16D4L20BE1P=FILL; ,= 

EDGE ; " 

Figure  A. 9   MAIN.BAS  (cont.) 
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1040  '  "Y" 

1050  X3=  X3  +  26 

1060  DRAW  "C=EDGE;BM=X3; ,=Y3;BU3  0BL1R5F9R4E9R5G14D16L4" 

1070  DRAW  "U16H14BR2BD1P=FILL; ,=EDGE;" 

1080  ■  "S" 

1090  X3=  X3  +  40 

1100  DRAW  "C=EDGE;BM=X3; ,=Y3;BU1U2E2R16H19U4E2R2  0F2D2G2" 

1110  DRAW  "L16F19D4G2L20H2BR2BU1P=FILL; ,=EDGE;" 

1120  '  "I" 

1130  X3=  X3  +  32 

1140  DRAW  "C=EDGE;BM=X3; ,=Y3 ;U30R4D30L4BE1P=FILL; , =EDGE ; " 

1150  '  "S" 

1160  X3=  X3  +  14 

117  0  DRAW  "C=EDGE;BM=X3; ,=Y3;BU1U2E2R16H19U4E2R2  0F2D2G2" 

1180  DRAW  "L16F19D4G2L20H2BR2BU1P=FILL; , =EDGE ; " 

1190  NEXT  I 

1200  LOCATE  22,4 

1210  INPUT  "WOULD  YOU  LIKE  INSTRUCTIONS  (Y/N)";Z$ 

1220  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  124  0 

1230  GOTO  2000 

1240  SCREEN  0 

1250  WIDTH  80 

12  60  COLOR  15,1 

127  0  CLS 

12  8  0  PRINT 

1290  PRINT  TAB (15) "THIS   IS  A  MENU   DRIVEN   PROGRAM  THAT 

WILL  ALLOW  YOU  TO" 
1300  PRINT  TAB (10) "EVALUATE  DIFFERENT   HELICOPTER  CONFIGU 

RATIONS   AND  DETERMINE" 
1310  PRINT  TAB (10) "THEIR  LIFT   AND  DRAG  COEFFICIENTS  AND 

THEIR  EQUIVALENT  FLAT" 

13  2  0  PRINT  TAB (10) "PLATE  AREA.   THE  DATA  GENERATED   DURING 

THE  TUNNEL  RUNS  WILL" 
1330  PRINT  TAB(10)"BE  USED  WITH  THIS  PROGRAM." 
1340  PRINT  TAB(15)"THE   FIRST   OPTION   IS   USED   TO  RECORD 

THE  DATA  DURING" 
1350  PRINT  TAB(10)"THE   TUNNEL  RUNS.   THE  SECOND   OPTION 

CONVERTS   THE  RECORDED" 
1360  PRINT  TAB(10)"RAW  COUNTS  TO   FORCES  AND  MOMENTS.   YOU 

WILL  HAVE  THE  CHANCE" 
1370  PRINT  TAB(10)"TO  PROVIDE  NAMES   FOR   EACH  OF  THESE 

FILES  AS  YOU  GO  ALONG." 
1380  PRINT  TAB(10)"THE  THIRD  OPTION   LETS  YOU  EXAMINE  THE 

FORCES  RECORDED.   IT" 
1390  PRINT  TAB(10)"WILL  BE  OF  NO  USE  TO  YOU  TO  EXAMINE  THE 

RAW  COUNTS . " 
1400  PRINT  TAB(15)"THE  FOURTH   OPTION  USES   THE  CONVERTED 

DATA  TO  COMPUTE" 
1410  PRINT  TAB(10)"THE   DESIRED   PARAMETERS.    THE  DATA 

MUST  BE  CONVERTED  USING" 

Figure  A. 9   MAIN.BAS  (cont.) 
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1420  PRINT  TAB (10) "OPTION  TWO  BEFORE  THE  PARAMETERS  CAN  BE 

CALCULATED.  YOU  WILL" 
1430  PRINT  TAB(10)"THEN  PROVIDE  A  NAME  FOR  THE   FILES  THAT 

STORE  THE  CALCULATED" 
1440  PRINT  TAB (10) "PARAMETERS.   USING   THESE  FILE   NAMES, 

THE  FIFTH  OPTION  LETS" 
1450  PRINT  TAB (10) "YOU  EXAMINE  THE  DATA  THAT  CAN  BE 

PLOTTED . " 

1455  PRINT: INPUT  "PRESS  RETURN  TO  CONTINUE" ;N 

1456  CLS 

1460  PRINT  TAB (15) "OPTION  SIX   MUST  BE  USED   TO  SORT  THE 

X  VALUES   BEFORE" 
1470  PRINT  TAB(10)"THEY  CAN  BE  PLOTTED.   THE  SEVENTH  OPTION 

LETS  YOU  ADD  A  GRID" 
1480  PRINT  TAB(10)"TO  THE   PLOT  AND  ALSO  MARK  THE  DATA 

POINTS.   THE  PLOTTING  IS" 
1490  PRINT  TAB(10)"IS   EXECUTED   WITH   THE   EIGHTH  OPTION. 

YOU  CAN  PLOT  UP  TO" 
1500  PRINT  TAB (10) "THREE  CURVES  PER  PLOT   AND  RECEIVE   A 

HARD  COPY  BY  PRESSING" 
1510  PRINT  TAB(10)"THE   CTRL-PRTSC   BUTTONS.    THE  FINAL 

OPTION  TERMINATES  THIS" 
152  0  PRINT  TAB (10) "PROGRAM." 
1530  PRINT:  INPUT  "PRESS  RETURN  TO  CONTINUE" ;N 

2000   ' PLOTTER  SET-UP 

2  010  SCREEN  0:  WIDTH  80:  COLOR  14,0,7:   CLS:   LOCATE  10,1 
202  0  PRINT  "WHICH  DISK  DRIVE  DO  YOU  WANT  TO  STORE  YOUR 

DATA  FILES  ON" ; 
2  03  0  PRINT 

2040  PRINT  "DRIVE  A  IS  EITHER  ON  YOUR  LEFT  OR  ON  TOP" 
2  050  PRINT  "DRIVE  B  IS  EITHER  ON  YOUR  RIGHT  OR  ON  THE 

BOTTOM" 
2  07  0  PRINT 

2  08  0  INPUT  "ENTER  [  A  OR  B  ]";FD$:  FD$=FD$+":"  :  PD$="C:" 
2090  GRD$="NO  GRID":  MRK$="MARK    " 
2100  OPT$=" CURVE" 
2110  GOSUB  2230 
2120  GOSUB  2400 

2130  X%=5+DX%:  Y%=12 :  LOCATE  X%,Y%,1,0,7 
2140  FIRST%=5:  LAST%=13 

2150  ANS$=INKEY$:  IF  ANS$=M"  THEN  2150 
2160  IF  ANS$=CHR$(0)+CHR$(80)  THEN  IF  X%<LAST%  THEN  X%= 

X%+1  ELSE  X%=FIRST% 
2170  IF  ANS$=CHR$(0)+CHR$(72)  THEN  IF  X%>FIRST%  THEN  X%= 

X%-1  ELSE  X%=LAST% 
2180  LOCATE  X%,Y%:   IF  ANS$<>CHR$ (13)  THEN  2150 
2190  FLAG1%=CSRLIN:   FLAGl%=FLAGl%-4 
2200  ON  FLAG1%  GOSUB  3120,2740,2630,2800,2680,3170,2870, 

3550,3680 
2210  IF  FLAG1%<>1  THEN  GOTO  2120  ELSE  GOTO  2130 

Figure  A. 9   MAIN.BAS  (cont.) 
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2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 

2420 
2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 


i INITIALIZING  CONSTANTS ■■ — — 

DIM  X(3  00),B(3  00),C(3  00),D(3  00),Y(300) 

DIM  FILESTK$(10) 

OUT  985,6 

SX%=5:   SY%=55 

M$="COMMAND  : " 

BLANK$="  " 

B$="  " 

SCREEN  0,1:  CLS 

HDPOS%=48 

KEY(l)  ON:  ON  KEY  (1)  GOSUB  2350 

RETURN 

i. CHANGE  FOREGROUND  COLORS 

FG%=FG%+1 

IF  FG%>15  THEN  FG%=1 

OUT  985, FG% 

RETURN 

i PRINT  MAIN  OPTION  MENU 

SCREEN  0,1:  COLOR  14,0,7:  CLS:  LOCATE  3,1 


PRINT  TAB (5 
********* 

PRINT 
PRINT  TAB (5 

PRINT  TAB (5 

PRINT  TAB (5 

DATA   *" 
PRINT  TAB (5 
*» 

PRINT  TAB (5 
NG      *" 
PRINT  TAB (5 
*•• 

PRINT  TAB (5 
*" 

PRINT  TAB (5 
*•• 

2510  PRINT  TAB (5 
*" 

252  0  PRINT 

2  53  0  PRINT  TAB (5 


"********        MAIN  OPTIONS  MENU 

"*  _  RECORD  TUNNNEL  DATA 

"*  REDUCE  RAW  TUNNEL  DATA 

"*  _  EXAMINE  DATA  FILE  FOR  TUNNEL 

"*  _  CALCULATE  Cd,  CI,  AND  E.F.P.A. 

"*  _  EXAMINE  DATA  FILE  FOR  PLOTTI 

"*  _  SORT  X  VALUES 

"*  _  PLOTTING  OPTION  MENU 

"*  DO  PLOTTING 

H-*  EXIT 


2540 
2550 
2560 
2570 
2580 
2590 
2600 


*************************************** 
********* 
LOCATE  3,65,1:  PRINT  "STATUS" 

LOCATE  4,55  :  PRINT  " " 

LOCATE  SX%,SY%  :  PRINT  "DATA  FILE  DRIVE    =  "+FD$ 


LOCATE  SX%+2,SY% 
LOCATE  SX%+3,SY% 
LOCATE  SX%+4,SY% 
LOCATE  SX%+5,SY% 


PRINT  "USER  OPTIONS 
PRINT  OPT$ 
PRINT  GRD$ 
PRINT  MRK$ 


Figure  A. 9   MAIN.BAS  (cont.) 
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2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
2710 
2720 
2730 
2740 
2750 
2760 
2770 
2780 
2790 
2800 
2810 
2820 
2830 
2840 
2850 
2860 
2870 
2880 

2890 
2900 

2910 

2920 

2930 

2940 

2950 
2960 

2970 
2980 
2990 
3000 

3010 


RETURN 

i DATA  EDITOR  FOR  TUNNEL  DATA— 

CHAINFILE$  =  PD$+"ADATA.BAS" 

DX%=2 

COMMON  FD$,PD$,DX% 

CHAIN  CHAINFILE$ 

RETURN 

i DATA  EDITOR  FOR  PLOTTING  DATA- 

CHAINFILE$  =  PD$+"BDATA.BAS" 

DX%=5 

COMMON  FD$,PD$,DX% 

CHAIN  CHAINFILE$ 

RETURN 

i DATA  REDUCTION 

CHAINFILE$  =  PD$+"RED.BASM 

DX%=3 

COMMON  PD$,FD$,DX% 

CHAIN  CHAINFILE$ 

RETURN 

i CALCULATE  ROUTINE 

CHAINFILE$  =  PD$  +  "COMP.BAS" 

DX%=4 

COMMON  PD$,FD$,DX% 

CHAIN  CHAINFILE$ 

REUTRN 

i PLOTTING  OPTION  MENU 

CLS:  SCREEN  0,1:  LOCATE  6,1 


PLOTTING  OPTION  MENU 

GRID  &     MARK 

GRID  &  NO  MARK 

NO  GRID  &     MARK 

NO  GRID  &  NO  MARK 

EXIT 


**  * 


PRINT  TAB(20)  "****** 

PRINT 

PRINT  TAB (20)  "* 
*» 

PRINT  TAB (2 00  "* 
*«• 

PRINT  TAB (20)  "* 
*" 

PRINT  TAB (20)  »* 
*» 

PRINT  TAB (20)  "* 

*H 

PRINT 

PRINT  TAB (20)  "************************************* 

X%=8:  Y%=27:  LOCATE  X%,Y%,1,0,7 

FIRST%=8:  LAST%=12 

ANS$=INKEY$:  IF  ANS$="n  THEN  2990 

IF  ANS$=CHR$(0)+CHR$(8  0)  THEN  IF  X%<LAST%  THEN  X%= 

X%+1  ELSE  X%=FIRST% 

IF  ANS$=CHR$(0)+CHR$(72)  THEN  IF  X%>FIRST%  THEN  X%= 

X%+1  ELSE  X%=LAST% 


Figure  A. 9   MAIN.BAS  (cont.) 
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3020 
3030 
3040 
3050 
3060 
3070 
3080 
3090 
3100 
3110 
3120 
3130 
3140 
3145 
3150 
3155 
3160 
3170 

3180 
3190 
3200 
3210 
3220 
3230 
3240 

3250 
3260 
3270 
3280 

3290 
3300 
3310 
3320 
3330 

3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 


LOCATE  X%,Y%:  IF  ANS$<>CHR$ (13)  THEN  2990 

FLAG%=CSRLIN :  FLAG%=FLAG%-7 

ON  FLAG%  GOSUB  3070,3080,3090,3100,3110 

DX%=7 

RETURN 

GRD$="GRID    ":  MRK$="MARK    "  :  RETURN 

GRD$="GRID    ":  MRK$="NO  MARK"  :  RETURN 

GRD$="NO  GRID":  MRK$="MARK    "  :  RETURN 

GRD$="NO  GRID":  MRK$="NO  MARK"  :  RETURN 

RETURN 

i RECORD  TUNNEL  DATA ■■ .— — 

CHAINFILE$  =  PD$  +  "RUNS.BAS" 

DX%=1 

COMMON  PD$,FD$,DX% 

CHAIN  CHAINFILE$ 

RETURN 

i SORTING 

CLS:  INPUT  "ENTER  NAME  OF  FILE  TO  BE  SORTED" ;FILE$ : 

FILE$=FD$+FILE$ 

PRINT  "READING  FILE  "+FILE$+" . . . " :  BEEP 

3370 

"SORTING" 

K%=NOD%-l:  R%=NOD% 

(L%<=R%) 

-1 

THEN  SWAP  X(J%) ,X(J%-1) :  SWAP 


THEN  SWAP  X(J%) ,X(J%-1) :  SWAP 


GOSUB 
PRINT 
L%=2: 
WHILE 

FOR  J%=R%  TO  L%  STEP 
IF  (X(J%-1)>X(J%) ) 
Y(J%) ,Y(J%-1) :  K%=J% 
NEXT 
L%=K%-1 

FOR  J%  =  L%  TO  R% 
IF  (X(J%-1)>X(J%) ) 
Y(J%) ,Y(J%-1) :  K%=J% 
NEXT 
R%=K%+1 
WEND 
CLS : BEEP 

PRINT  "SORTED  FILE  ";FILE$;"  BEING  SAVED 
3460 
DX%=6 
RETURN 

i READING  A  FILE 

OPEN  FILE$  FOR  INPUT  AS  #2 

INPUT  #2 , NOD% , Yl$ , Xl$ , CONFIG$ , CONF 

INPUT  #2,Q 

FOR  J  =  1  TO  NOD% 

INPUT  #2,Y(J) ,X(J) 

NEXT 

#2 


GOSUB 


CLOSE 
RETURN 


■SAVE  A  FILE- 


Figure  A. 9   MAIN.BAS  (cont.) 
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34  60  OPEN  FILE$  FOR  OUTPUT  AS  #1 

3470  WRITE  #l,NOD%, Yl$ ,X1$ , CONFIG$ , CONF 

3480  WRITE  #1,Q 

3490    FOR  J  =  1  TO  NOD% 

3500        WRITE  #1,Y(J),X(J) 

3510        NEXT 

3520  CLOSE  #1 

3530  RETURN 

3540  ' INITIALIZING  PLOTTER 

3  550  CLS:  INPUT  "NUMBER  OF  DATA  FILE(S)  TO  BE  PLOTTED  ON 

THE  SAME  PLOT  =";NOF% 
3  560  PRINT  "ENTER  FILE  NAME  (S)  : " : BEEP 
3570  FOR  1%  =  1  TO  NOF% 
3580    PRINT  "FILE  #";: PRINT  USING  "##";I%;  :  INPUT  "  =" ; 

FILE$ 
3590    FILE$  =  FD$  +  FILE$ 
3600    FILESTK$(I%)  =  FILE$ 
3  601    PRINT  "HOW  DO  YOU  WANT  THIS  CURVE  PLOTTED  ?   ENTER 

(POINT)  TO  JUST  PLOT" 
3602    INPUT  "THE  POINTS  OR  (CURVE)  TO  PRODUCE  A  CURVE 

FIT";OPT$ 
3610    CURVE$(I%)  =  OPT$ 
3  620    NEXT 

3630  PFILE$  =  PD$+"MAIN.BAS":  FILE2$=  PD$  +  "PLOTTER. BAS" 
3640  SL%  =  100 
3  650  CHAIN  FILE2$ , 30 , ALL 
3660  DX%=8 
3  670  RETURN 
3680  CLS: 
3  690  LOCATE  10,5:  PRINT  "PROGRAM  TERMINATED,  REMOVE  DATA 

FILES  FROM  DISK  DRIVE" 
37  00  LOCATE  11,5:  PRINT  "AND  SECURE  COPMUTER  AND  AMPLIFIE 

S" 
3710  END 

Figure  A. 9   MAIN. BAS  (cont.) 
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NAME:   Data  Acquisition  And  Control  (DAAC) 
HEADER  for  BASICA 

FILE  NAME:   DACHDR.BAS 

DOS  DEVICE  NAME:   DAAC 

RESERVED  FUNCTION  NAMES: 

AINM,  AINS,  AINSC,  AOUM,  AOUS , 
BINM,  BINS,  BITINS,  BITOUS ,  BOUM,  BOUS , 
CINM,  CINS,  CSET,  DELAY 
RESERVED  DEF  SEG  VALUE  NAME:   DSEG 

NAMES  DEFINED  AND  USED  BY  HEADER: 

ADAPT%,  AI,  COUNT,  FOUND% , 
HNAME$,  SG%,  STAT% 


When  using  the  BASICA  Interpreter,  this  header 
must  be  executed  before  any  function  calls  are 
made  that  access  the  DAAC  adapter.   It  initializes 
a  number  of  variables  for  each  function  call.   These 
variables  are  reserved  and  should  not  be  used  except 
to  access  the  DAAC  adapter.   This  routine  also  does  a 
DEF  SEG  to  the  segment  where  the  DAAC  Device  Driver 
(DAC.COM)  is  loaded.   If  you  execute  a  DEF  SEG  to 
access  other  hardware,  you  must  DEF  SEG  to  the  segment 
of  the  DAAC  Device  Driver  before  any  subsequent 
calls  to  access  the  DAAC  adapter. 


FOUND%  =  0 
SG%  =  &H2E 

' Start  searching  the  interrupt  vectors  until  you  find 
'one  that  points  to  the  DAAC  device  driver. 
'Do  a  DEF  SEG  to  that  segment. 
WHILE  ((SG%  <=  &H3E)  AND  (FOUND%  =  0)) 

DEF  SEG  =  0 

DSEG  =  PEEK(SG%)  +  PEEK(SG%  +  1)  *  256 

DEF  SEG  =  DSEG 

HNAME$=',M 

FOR  AI=10  TO  17 

HNAME$  =  HNAME$  +  CHR$ (PEEK (AI) ) 

NEXT  AI 

IF  HNAME$  =  "DAAC     "  AND  PEEK (18)  +  PEEK (19)  <> 

0  THEN  F0UND%  =  1 

SG%  =  SG%  +  4 


WEND 

IF  FOUND%  = 

NOT  FOUND"  : 


0  THEN  PRINT  "ERROR: 
END 


DEVICE  DRIVER  DAC.COM 


Figure  A. 10   RUNS . BAS  -  Data  Acquisition  Program 
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•Now  initialize  all  function  name 
'to  access  the  device  driver. 


variables  for  calls 


AINM         =  PEEK(&H13)  *  256  + 

AINS         =  PEEK(&H15)  *  256  + 

AINSC        =  PEEK(&H17)  *  256  + 

AOUM         =  PEEK(&H19)  *  256  + 

AOUS         =  PEEK(&H1B)  *  256  + 

BINM         =  PEEK(&H1D)  *  256  + 

BINS         =  PEEK(&H1F)  *  256  + 

BITINS       =  PEEK(&H21)  *  256  + 

BITOUS       =  PEEK(&H23)  *  256  + 

BOUM         =  PEEK(&H25)  *  256  + 

BOUS         =  PEEK(&H27)  *  256  + 

CINM  •  =  PEEK(&H29)  *  256  + 

CINS         =  PEEK(&H2B)  *  256  + 

CSET         =  PEEK(&H2D)  *  256  + 

DELAY        =  PEEK(&H2F)  *  256  + 

•Finally,  execute  any  call  to 
•device  driver  from  any  former  invocation  of  BASIC. 
ADAPT%  =  0 

COUNT  =  1  . 
STAT%  =  0 
CALL  DELAY  ( ADAPT %,  COUNT,  STAT%) 


PEEK(&H12) 
PEEK(&H14) 
PEEK(&H16) 
PEEK(&H18) 
PEEK(&H1A) 
PEEK(&H1C) 
PEEK(SHIE) 
PEEK(&H2  0) 
PEEK(&H22) 
PEEK(&H24) 
PEEK(&H2  6) 
PEEK(&H28) 
PEEK(&H2A) 
PEEK(&H2C) 
PEEK(&H2E) 
re-initialize  the 


'End  of  DAAC  BASICA  Header 
i 

REM —  RUNS.BAS;   PROGRAM  TO  RECORD  THE  DATA  FROM 
WIND 

REM —  TUNNEL  RUNS 
KEY  OFF:  COLOR  15,1,4:  CLS 
CLS 

PRINT:  PRINT 

PRINT  TAB (10)  "HELO  CONFIGURATIONS" 
PRINT:  PRINT 


PRINT  TAB (5) 
PRINT  TAB (5) 
PRINT  TAB (5) 
PRINT  TAB (5) 


PRINT 
PRINT 
PRINT 
PRINT 


TAB ( 5 ) 
TAB ( 5 ) 
TAB ( 5 ) 
TAB  (5 


1.  ATTACK  NOSE, 

2.  ATTACK  NOSE, 

3.  ATTACK  NOSE, 

4.  SMOOTH  NOSE, 

5.  SMOOTH  NOSE, 

6.  SMOOTH  NOSE, 

7.  BLUNT  NOSE, 

8.  BLUNT  NOSE, 

9.  BLUNT  NOSE, 


STRAIGHT  TAIL" 

LOW  TAIL" 

HIGH  TAIL" 

STRAIGHT  TAIL" 

LOW  TAIL" 

HIGH  TAIL" 
STRAIGHT  TAIL" 
LOW  TAIL" 
HIGH  TAIL" 


PRINT  TAB (5) 

PRINT :  PRINT 

INPUT  "WHICH  CONFIGURATION  IS  BEING  RUN";N 

PRINT 

INPUT  "DOES  CONFIGURATION  INCLUDE  LANDING  GEAR";Z$ 

ON  N  GOSUB  4550,4620,4690,4750,4820,4870,4920,4990, 


5040 
i 


WEIGHTS  ARE  IN  POUNDS' 


Figure  A. 10   RUNS.BAS  (cont.) 
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960    IF  CONF  =  10   THEN  W  =  22.7892 
970    IF  CONF  =  15   THEN  W  =  23.4623 
980    IF  CONF  =  20   THEN  W  =  24.6945 
990    IF  CONF  =  25   THEN  W  =  25.3676 
1000    IF  CONF  =  30   THEN  W  =  24.6945 
1010    IF  CONF  =  35   THEN  W  =  25.3676 
1020    IF  CONF  =  40   THEN  W  =  20.8134 
1030    IF  CONF  =  45   THEN  W  =  21.5866 
1040    IF  CONF  =  50   THEN  W  =  22.7187 
1050    IF  CONF  =  55   THEN  W  =  23.4913 
1060    IF  CONF  =  60   THEN  W  =  22.7187 
1070    IF  CONF  =  65   THEN  W  =  23.4913 
1080    IF  CONF  =  70   THEN  W  =  21.3315 
1090   '  IF  CONF  =  75   THEN  W  =  23.2892 
1100    IF  CONF  =  80   THEN  W  =  23.2368 
1110    IF  CONF  =  85   THEN  W  =  25.1945 
1120    IF  CONF  =  90   THEN  W  =  23.2368 
1130    IF  CONF  =  95   THEN  W  =  25.1945 
1140    CLS 

1150    REM RECORD  TUNNEL  DATA — ■ 

1160  DIM  L(100) ,D(100) ,Y(100) ,PM(100) ,YM(100) ,RM(100) ,AOA( 

100) 
1170  DIM  DAT(399) ,DAT%(399) ,DAT1(399) ,DAT1%(399) 
118  0  COLOR  15,1:  KEY  OFF:  CLS 

1190  PRINT  "RECORDING  OF  WIND  TUNNEL  RAW  DATA" 
1200  INPUT  "WHAT  IS  THE  TUNNEL  SPEED  (Q)  FOR  THIS  RUN";Q' 
1210  PRINT 
2  67  0  PRINT  "WITHOUT  THE  MODEL  IN  PLACE,  ADJUST  THE  ZERO 

SET  SCREWS  TO  ZERO" 
2  671  PRINT  "OUT  EACH  AMPLIFIER.   AFTER  ALL  AMPLIFIERS  ARE 

ZEROED,  LOAD  THE" 
2  67  2  PRINT  "MODEL  ONTO  THE  STING.   WHEN  THE  MODEL  IS 

MOUNTED,  PRESS  RETURN" 
2673  INPUT  "TO  RECORD  A  NO  FORCE  ZERO  READING" ;N 
2  690  PRINT: PRINT  "   ZD        ZL       ZY         ZPM       ZYM 

ZRM      ZAOA" 
2700  STAT%=0:   MODE%=0:  STOR%=0 :  COUNT=100:  RATE=500 
2710  ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=3 
272  0  CALL  AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% ,MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 
2730  ZD  =  0:ZPM=0:ZL=0:ZYM=0 
2740  IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT% : END 
2750  FOR  J  =  0  TO  396  STEP  4 
2760  DAT(J)=(DAT%(J)/204.8)-10 
2770  ZD  =ZD  +  DAT(J) 
2780  NEXT  J 
2790  ZD  =  ZD/100 
2800  FOR  J  =  1  TO  397  STEP  4 
2810  DAT(J)=(DAT%(J)/204.8)-10 

Figure  A. 10   RUNS . BAS  (cont.) 
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ZL  =  ZL  +  DAT (J) 

NEXT  J 

FOR  J  =  2  TO  398  STEP  4 

DAT (J) = (DAT% (J) /204 . 8 ) -10 

ZPM  =  ZPM  +  DAT (J) 

NEXT  J 

FOR  J  =  3  TO  399  STEP  4 

DAT(J)=(DAT%(J)/204.8)-10 

ZYM  =  ZYM  +  DAT (J) 

NEXT  J 

STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=2 

CALL  AINSC ( ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 )  , STAT% ) 

ZY  =0 : ZRM-0 : ZAOA=0 

IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT% : END 

FOR  J  =  0  TO  297  STEP  3 

DAT1 (J) = (DAT1% ( J) /2  04 . 8 ) -10 

ZRM  =  ZRM  +  DAT1(J) 

NEXT  J 

ZRM  =  ZRM/100 

FOR  J  =  1  TO  298  STEP  3 

DAT1 (J) = (DAT1% ( J) /2  04 . 8 ) -10 

ZY  =  ZY  +  DAT1(J) 

NEXT  J 

FOR  J  =  2  TO  299  STEP  3 

DATl(J)=(DATl%(J)/2  04.8)-10 

ZAOA  =  ZAOA  +  DAT1(J) 

NEXT  J 

ZL=ZL/100 : ZPM=ZPM/100 : ZYM=ZYM/100 : ZY=ZY/100 : ZAOA=ZAOA 

/100 

LOCATE  11,1:  PRINT  USING  »+# . ###" ;ZD:  LOCATE  11,10: 

PRINT  USING  "+#.###" ;ZL 

LOCATE  11,19:  PRINT  USING  "+#.###" ;ZY 

LOCATE  11,28:  PRINT  USING  "+#.###"; ZPM 

LOCATE  11,37:  PRINT  USING  "+#.###"; ZYM 

LOCATE  11,46:  PRINT  USING  "+#.###"; ZRM 

LOCATE  11,55:  PRINT  USING  "+#.###"; ZAOA 

REM 

PRINT: PRINT 

INPUT  "AFTER  PLACING  ALL  CAL  SWITCHES  TO  +  SETTING  HIT 

RETURN" ;X 

PRINT 

PRINT  "  CALD      CLL     CALY      CALPM      CALYM 

CALRM     CALAOA" 

STAT%=0:   MODE%=0:  STOR%=0 :  COUNT=100:  RATE=500 

ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=3 

CALL  AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 


Figure  A. 10   RUNS . BAS  (cont.) 
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CALD=  0 : CLL=  0 : CALYM=  0 : CALPM=  0 

IF  STAT%<>  0,THEN  PRINT  USING  "EXECUTION  ERROR  ###" 
; STAT% : END 

FOR  J  =  0  TO  396  STEP  4 
DAT (J) = (DAT% (J) /2  04 . 8 ) -10 
CALD  =CALD  +  DAT (J) 
NEXT  J 

CALD  =  CALD/ 100 
FOR  J  =  1  TO  397  STEP  4 
DAT (J) = (DAT% (J) /204 . 8 ) -10 
CLL  =  CLL  +  DAT (J) 
NEXT  J 

FOR  J  =  2  TO  398  STEP  4 
DAT (J) = (DAT% (J) /2  04 . 8 ) -10 
CALPM  =  CALPM  +  DAT (J) 
NEXT  J 

FOR  J  =  3  TO  399  STEP  4 
DAT(J)=(DAT%(J)/2  04.8)-10 
CALYM  =  CALYM  +  DAT (J) 
NEXT  J 

STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 
ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=2 
CALL  AINSC ( ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 
CALY  =0 : CALRM=0 : CALAOA=0 

IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 
STAT%:END 

FOR  J  =  0  TO  297  STEP  3 
DAT1 (J) = (DAT1% ( J) /2  04 . 8) -10 
CALRM  =  CALRM  +  DAT1(J) 
NEXT  J 

CALRM  =  CALRM/ 100 
FOR  J  =  1  TO  298  STEP  3 
DAT1 (J)  =  (DAT1% (J) /2  04 . 8 ) -10 
CALY  =  CALY  +  DAT1(J) 
NEXT  J 

FOR  J  =  2  TO  299  STEP  3 
DAT1 ( J) = (DAT1% ( J) /2  04 . 8 ) -10 
CALAOA  =  CALAOA  +  DAT1(J) 
NEXT  J 

CLL=CLL/100 : CALPM=CALPM/ 1 0 0 : CALYM=CALYM/100 : CALY=CALY 
/ 1 0  0 : CALAOA=CALAOA/ 100 

LOCATE  18,1:  PRINT  USING  "+#.###"; CALD 
LOCATE  18,10:  PRINT  USING  "+# 
19:  PRINT  USING  "  +  # 

"  +  # 

"  +  # 

"  +  # 

"  +  # 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
INPUT 


18, 
18, 
18, 
18, 
18, 


28 
37 

46 
55 


"PRESS 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
RETURN 


USING 
USING 
USING 
USING 
USING 
USING 


###" 
###" 
###" 
###" 
###" 
###" 


CLL 

CALY 

CALPM 

CALYM 

CALRM 

CALAOA 


TO  CONTINUE" ;N 


Figure  A. 10   RUNS . BAS  (cont. ) 
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CLS 

PRINT:  PRINT  "REPLACE  THE  CAL  SWITCHES  TO  THE  CENTER 

POSITION" 

PRINT  "AFTER  YOU  HAVE  GOTTEN  THE  WIND  TUNNEL  UP  TO 

SPEED  AND  ARE  READY" 

PRINT  "TO  RECORD  DATA  PRESS  RETURN.   ONCE  THE  HEADINGS 

ARE  PRINTED  THE" 

PRINT  "F2  KEY  WILL  RECORD  THE  DATA.   THE  Fl  KEY  WILL 

SAVE  THE  DATA  AND  " 

PRINT  "RERUN  YOU  TO  THE  MAIN  MENU.   YOU  CAN  TAKE  READ 

INGS  FOR  ANY  ANGLE" 

PRINT  "OF  ATTACK  BETWEEN  +10  AND  -8  DEGREES." 

PRINT 

PRINT  "THIS  PROGRAM  CORRECTS  FOR  DRIFT  IN  THE  BALANCE 

AND  AMPLIFIERS." 

PRINT  "YOUR  FIRST  AND  LAST  DATA  POINTS  SHOULD  BE  TAKEN 

AT  ZERO  ANGLE  OF" 

PRINT  "ATTACK.   ALSO,  DO  NOT  TAKE  MORE  THAN  ONE  DATA 

POINT  FOR  THE  SAME  " 

PRINT  "ANGLE  OF  ATTACK  EXCEPT  FOR  YOUR  FIRST  AND  LAST 

DATA  POINT.   YOU" 

PRINT  "ARE  LIMITED  TO  100  .DATA  POINTS.   REMEMBER,  EACH 

RUN  IS  FOR  ONE" 

PRINT  "SPECIFIC  'Q'  SETTING" 

INPUT  "PRESS  RETURN  TO  CONTINUE" ;N 

CLS 

PRINT 

ROLL 

PRINT  " 

MOM. 
SOAP  =  0: 
FOR  K  =  1 
ON  KEY(l) 
ON  KEY (2) 
KEY(l)  ON 
IF  SOAP  = 
IF  SOAP  = 
GOTO  3720 
SOAP  =  0 
NOD%  =  K 
NEXT  K 
GOTO  4370 

REM   SET  STOP  FLAG 
SOAP  =  1 
RETURN 

REM  STEPS  TO  RECORD  DATA 

STAT%=0:   MODE%=0:  STOR%=0 :  COUNT=100:  RATE=500 
ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=3 
CALL  AINSC ( ADAPT % , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 
, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 


DRAG 
AOA 


LIFT      YAW 
TUNNEL" 

SPEED" 


PITCH 


MOM. 


YAW 


MOM. 


N=4 

TO  100 
GOSUB  38  2  0 
GOSUB  3  850 
KEY ( 2 )  ON 
2  THEN  GOTO  3780 
1  THEN  GOTO  4  370 


'SET  STOP  FLAG 
'RECORD  DATA 


Figure  A. 10   RUNS . BAS  (cont.) 
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3890  D(K)=0:L(K)=0:YM(K)=0:PM(K)=0 

3900  IF  STAT%<>  0,THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT% : END 
3910  FOR  J  =  0  TO  396  STEP  4 
3920  DAT(J)=(DAT%(J)/204.8)-10 
3930  D(K)  =D(K)  +  DAT(J) 
3940  NEXT  J 
3950  D(K)  =  D(K)/100 
3960  FOR  J  =  1  TO  397  STEP  4 
3970  DAT(J)=(DAT%(J)/2Q4.8)-10 
3980  L(K)  =  L(K)  +  DAT(J) 
3990  NEXT  J 

4000  FOR  J  =  2  TO  398  STEP  4 
4010  DAT(J)=(DAT%(J)/204.8)-10 
4020  PM(K)  =  PM(K)  +  DAT(J) 
403  0  NEXT  J 

4040  FOR  J  =  3  TO  399  STEP  4 

4041  DAT(J)=(DAT%(J)/204.8)-10 
4050  YM(K)  =  YM(K)  +  DAT(J) 
4060  NEXT  J 

4070  STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

4080  ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=2 

4  090  CALL  AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT%) 
4100  Y(K)  =0:RM(K)=0:  AOA(K)=0 
4110  "IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT% : END 
4120  FOR  J  =  0  TO  297  STEP  3 
4130  DATl(J)=(DATl%(J)/204.8)-10 
4140  RM(K)  =  RM(K)  +  DAT1(J) 
4150  NEXT  J 

4160  RM(K)  =  RM(K)/100 
4170  FOR  J  =  1  TO  298  STEP  3 
4180  DATl(J)=(DATl%(J)/204.8)-10 
4190  Y(K)  =  Y(K)  +  DAT1(J) 
4200  NEXT  J 

4210  FOR  J  =  2  TO  299  STEP  3 
4220  DATl(J)=(DATl%(J)/204.8)-10 
4230  AOA(K)  =  AOA(K)  +  DAT1(J) 
4240  NEXT  J 
4250  AOA(K)=AOA(K)/100:L(K)=L(K)/100:PM(K)=PM(K)/100:YM(K) 

=YM(K)/100:Y(K)=Y(K)/100 

4255  B  =  AOA(K)  -  ZAOA 

4256  AOA(K)  =  (5.8469*B)  +  (. 0077583* (BA2) ) 

4  2  60  IF  N>  2  3  THEN  N=4 : CLS : PRINT  "  DRAG      LIFT       YAW 
PITCH      YAW      ROLL      AOA         TUNNEL" : PRINT  " 
MOM.       MOM.     MOM.  SPEED" 

4270  LOCATE  N,l:  PRINT  USING  "+#.###" ;D (K) 
4280  LOCATE  N,10:  PRINT  USING  "+#.###" ;L(K) 
4290  LOCATE  N,19:  PRINT  USING  "+#.###" ;Y (K) 
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4300  LOCATE  N,28:  PRINT  USING  "  +  # . ###" ;PM(K) 
4310  LOCATE  N,37:  PRINT  USING  "+#.###" ;YM(K) 
4320  LOCATE  N,46:  PRINT  USING  "+#.###" ;RM(K) 
4330  LOCATE  N,55:  PRINT  USING  »+#.###" ;AOA(K) 
43  4  0  N=N+1 
4350  SOAP  =  2 

43  60  RETURN 

4370  SOAP  =  0 

4371  ■ CORRECT  FOR  DRIFT 

4372  DIFFL  = (L(NOD%) -L(l) ) :  DIFFD= (D (NOD%) -D(l) ) : DIFFY= (Y ( 
NOD%)-Y(l) ) 

4373  DIFFPM=(PM(NOD%) -PM(1) ) : DIFFYM= (YM(NOD%) -YM(1) ) :DIFFR 
M=(RM(NOD%)-RM(l) ) : DIFFAOA= (AOA(NOD%) -AOA(l) ) 

4374  D=NOD%-l 

4375  CORL=DIFFL/D:CORD=DIFFD/D:CORY=DIFFY/D:CORPM=DIFFPM/D 
: CORYM=DIFFYM/D : CORRM=DIFFRM/D : CORAOA=DIFFAOA/D 

4376  A=l 

4377  FOR  K  =  2  TO  NOD% 

4378  L(K)=L(K)-(A*CORL) 

4379  D(K)=D(K)-(A*CORD) 

4380  Y(K)=Y(K)-(A*CORY) 

4381  PM(K)=PM(K)-(A*CORPM) 

4382  YM(K)=YM(K)-(A*CORYM) 

4383  RM(K)=RM(K)-(A*CORRM) 

4385  A  =  A+l 

4386  NEXT  K 

4387  ■ SAVE  DATA  TO  FILE 

4  390  CLS:  PRINT 

4400  INPUT  "WHAT  IS  THE  NAME  FOR  THE  FILE  TO  STORE  THE  RAW 

DATA";FILE$ 
4410  FILE$=FD$+FILE$ 

44  2  0  OPEN  FILE$  FOR  OUTPUT  AS  #1 
4421  WRITE  #l,NOD%,CONFIG$,CONF 
4430  WRITE  #1,Q,W 

4470  WRITE  #1 , ZL, ZD, ZY , ZPM, ZYM, ZRM, ZAOA 

4480  WRITE  #1 , CLL, CALD, CALY, CALPM, CALYM, CALRM, CALAOA 

4490     FOR  J=l  TO  NOD% 

4500       WRITE  #lfL(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,AOA(J) 

4510       NEXT  J 

4520  CLOSE  #1 

4521  CLS: PRINT  "DO  YOU  WANT  TO  TAKE  ANOTHER  SET  OF  DATA 
POINTS  FOR  THE  SAME" 

4  52  2  INPUT  "CONFIGURATION  BUT  A  DIFFERENT  'Q'  SETTING  (Y/N 
) ";ANS$ 

4523  IF  ANS$="Y"  OR  ANS$="y"  THEN  GOTO  4524  ELSE  GOTO  4530 

4524  PRINT:  INPUT  "WHAT  IS  THE  'Q'  VALUE  FOR  THIS  RUN";Q 

4525  GOTO  3610 

453  0  ■ RETURN  TO  MAIN  PROGRAM 

4531  COMMON  PD$ , FD$ , DX% 

453  2  CHAIN  PD$+"MAIN. BAS" , 2110 
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454  0  REM SUBROUTINES  FOR  CONFIGURATION  MARKING— 

4550  IF  Z$  =  "Y"  OR  Z$  -  "y"  THEN  GOTO  4590 

4560  CONF  =  10 

4  570  CONFIG$=" ATTACK  NOSE/STRAIGHT  TAIL  WITHOUT  GEAR" 

4580  GOTO  4610 

4590  CONF  =  15 

4600  CINFIG$=" ATTACK  NOSE/STRAIGHT  TAIL  WITH  GEAR" 

4610  RETURN 

4620  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4660 

4630  CONF  =  20 

4640  CONFIG$=" ATTACK  NOSE/LOW  TAIL  WITHOUT  GEAR" 

4650  GOTO  4680 

4660  CONF  =  25 

4670  CINFIG$=" ATTACK  NOSE/LOW  TAIL  WITH  GEAR" 

4  680  RETURN 

4690  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4730 

4700  CONF  =  30 

4710  CONFIG$=" ATTACK  NOSE/HIGH  TAIL  WITHOUT  GEAR" 

4720  GOTO  4740 

4730  CONF  =  35:  CINFIG$="ATTACK  NOSE/HIGH  TAIL  WITH  GEAR" 

4740  RETURN 

4750  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4790 

4760  CONF  =  40 

4770  CONFIG$="SMOOTH  NOSE/STRAIGHT  TAIL  WITHOUT  GEAR" 

4780  GOTO  4810 

4790  CONF  =  45 

4800  CINFIG$="SMOOTH  NOSE/STRAIGHT  TAIL  WITH  GEAR" 

4810  RETURN 

4820  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4850 

4830  CONF=50:CONFIG$=" SMOOTH  NOSE/LOW  TAIL  WITHOUT  GEAR" 

4840  GOTO  4860 

4850  CONF  =  55:  CINFIG$="SMOOTH  NOSE/LOW  TAIL  WITH  GEAR" 

4860  RETURN 

4870  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4900 

4880  CONF=60:CONFIG$="SMOOTH  NOSE/HIGH  TAIL  WITHOUT  GEAR" 

4890  GOTO  4910 

4900  CONF  =  65:  CINFIG$="SMOOTH  NOSE/HIGH  TAIL  WITH  GEAR" 

4910  RETURN 

4920  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  4960 

4930  CONF  =  70 

4940  CONFIG$="ATTACK  NOSE/STRAIGHT  TAIL  WITHOUT  GEAR" 

4950  GOTO  4980 

4960  CONF  =  75 

497  0  CINFIG$=" ATTACK  NOSE/STRAIGHT  TAIL  WITH  GEAR" 

498  0  RETURN 

4990  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  5020 

5000  CONF  =  80:CONFIG$="BLUNT  NOSE/LOW  TAIL  WITHOUT  GEAR" 

5010  GOTO  5030 

5020  CONF  =  85:  CINFIG$="BLUNT  NOSE/LOW  TAIL  WITH  GEAR" 

503  0  RETURN 
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5040  IF  Z$  =  "Y"  OR  Z$  =  "y"  THEN  GOTO  5070 

5050  CONF=90:CONFIG$=" BLUNT  NOSE/HIGH  TAIL  WITHOUT  GEAR" 

5060  GOTO  5080 

5070  CONF  =  95:  CINFIG$="SMOOT  NOSE/HIGH  TAIL  WITH  GEAR" 

5080  RETURN 
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10    REM RED.BAS     (CONVERT  RAW  COUNTS  TO  FORCES) 

20  REM 
3  0  REM 
40    COLOR  15,1: KEY  OFF:  CLS 

45  DIM  L(100)  ,L1(100)  ,1*2(100)  ,D(100)  ,D1(100)  ,D2(100)  , 

(100) ,Y1(100) ,Y2(100) 

46  DIM  PM(IOO) ,PM1(100) ,PM2(100) ,YM(100) ,YM1(100) ,YM2(10 
0) ,RM(100) ,RM1(100) ,RM2(100) 

47  DIM  LF(IOO) ,DF(100) ,AOA(100) 

50    PRINT: PRINT "WHAT  IS  THE  NAME  OF  THE  FILE  THAT  CONTAINS 

THE  RAW  DATA" 
60    INPUT"THAT  YOU  WISH  TO  CONVERT" ; FILE $ 
70    FILE$=FD$+FILE$ 

8  0    • READ  IN  RAW  DATA  FILE — 

9  0  OPEN  FILE$  FOR  INPUT  AS  #1 
100  INPUT  #l,NOD%,CONFIG$,CONF 
110   INPUT  #1,Q,W 

150   INPUT  #l,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 

160   INPUT  #l,CLL,CALD,CALY,CALPM,CALYM,CALRM,CALAOA 

170      FOR  J  =  1  TO  NOD% 

180      INPUT  #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,AOA(J) 

190      NEXT  J 

2  00   CLOSE  #1 

210  PRINT: PRINT "PLEASE  WAIT  WHILE  PERFORMING  CONVERSION" 

211  PRINT  "THE  CONVERSION  WILL  TAKE  APPROXIMATELY  TEN  MIN 
UTES  FOR  " 

212  PRINT  "TWENTY  DATA  POINTS" 

220  ' READ  IN  CALIBRATION  CONSTANTS 

221  OPEN  "C: CONST"  FOR  INPUT  AS  #1 

222  INPUT  # 1 , INCALL , INCALD , INCALY , INCALPM , INCALYM , INCALRM 
240   INPUT  #1,K1LP0S,K2LP0S,K1DP0S,K2DP0S,K1YP0S,K2YP0S 

2  50   INPUT  #1,K1PMP0S,K2PMP0S,K1YMP0S,K2YMP0S,K1RMP0S, 
K2RMPOS 

2  60   INPUT  #1,K1LNEG,K2LNEG,K1DNEG,K2DNEG,K1YNEG,K2YNEG 
27  0   INPUT  #1,K1PMNEG,K2PMNEG,K1YMNEG,K2YMNEG,K1RMNEG, 

K2RMNEG 
280   INPUT  #1,DDDL1P,DDDL2P,DYDL1P,DYDL2P,DPMDL1P,DPMDL2P, 

DYMDL1P,DYMDL2P 
290   INPUT  #1,DRMDL1P,DRMDL2P 
300   INPUT  #1,DLDD1P,DLDD2P,DYDD1P,DYDD2P,DPMDD1P,DPMDD2P, 

DYMDD1P,DYMDD2P 
310   INPUT  #1,DRMDD1P,DRMDD2P 
320   INPUT  #1,DLDY1P,DLDY2P,DDDY1P,DDDY2P,DPMDY1P,DPMDY2P, 

DYMDY1P,DYMDY2P 
330   INPUT  #1,DRMDY1P,DRMDY2P 
340   INPUT  #1,DLDPM1P,DLDPM2P,DDDPM1P,DDDPM2P,DYDPM1P,DYDP 

M2P,DYMDPM1P,DYMDPM2P 

3  50   INPUT  #1,DRMDPM1P,DRMDPM2P 

360   INPUT  #1,DLDYM1P,DLDYM2P,DDDYM1P,DDDYM2P,DYDYM1P,DYDY 
M2  P , DPMDYM1P , DPMDYM2  P 
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INPUT  #1,DRMDYM1P,DRMDYM2P 

INPUT  #1, DLDRM1P, DLDRM2P, DDDRM1P, DDDRM2P, DYDRM1P, DYDR 

M2  P , DPMDRM1P , DPMDRM2  P 

INPUT  # 1 , DYMDRM1P , DYMDRM2  P 

INPUT  # 1 , DDDL1N , DDDL2N , DYDL1N , DYDL2N , DPMDL1N , DPMDL2N , 

DYMDL1N,DYMDL2N 

INPUT  #1,DRMDL1N,DRMDL2N 

INPUT  # 1 , DLDD1N , DLDD2N , DYDD1N , DYDD2N , DPMDD1N , DPMDD2N , 

DYMDD1N,DYMDD2N 

INPUT  #1,DRMDD1N,DRMDD2N 

INPUT  # 1 , DLDY1N , DLDY2N , DDDY1N , DDDY2N , DPMDY1N , DPMDY2N , 

DYMDY1N,DYMDY2N 

INPUT  #1,DRMDY1N,DRMDY2N 

INPUT  # 1 , DLDPM1N , DLDPM2N , DDDPM1N , DDDPM2N , DYDPM1N , DYDP 

M2N,DYMDPM1N 

INPUT  # 1 , DYMDPM2N , DRMDPM1N , DRMDPM2N 

INPUT  # 1 , DLDYM1N , DLDYM2N , DDDYM1N , DDDYM2N , DYDYM1N , DYDY 

M2N,DPMDYM1N 

INPUT  #1,DPMDYM2N,DRMDYM1N,DRMDYM2N 

INPUT  # 1 , DLDRM1N , DLDRM2N , DDDRM1N , DDDRM2N , DYDRM1N , DYDR 

M2N,DPMDRM1N 

INPUT  # 1 , DPMDRM2N , DYMDRM1N , DYMDRM2N 

CLOSE  #1 

1 CONVERT  RAW  COUNTS  TO  FORCES 

FOR  J  =  1  TO  NOD% 

A=  ( (INCALL/(CLL-ZL) )*(L(J)-ZL) ) 
B  =  ( (INCALD/(CALD-ZD) ) *(D(J)-ZD) ) 
C  =  ( (INCALY/(CALY-ZY) )*(Y(J)-ZY) ) 
D=( (INCALPM/(CALPM-ZPM) ) * (PM(J) -ZPM) ) 
E=( (INCALYM/(CALYM-ZYM) ) * (YM(J) -ZYM) ) 
F=( (INCALRM/(CALRM-ZRM) ) * (RM( J) -ZRM) ) 

IF  A  <=  0  THEN  K1=K1LNEG:K2=K2LNEG 
:K2=K2LPOS 

L(J)  =  (K1*A)  +  (K2*AA2) 

IF  B  <=  0  THEN  K1=K1DNEG:K2=K2DNEG 
:K2=K2DPOS 

D(J)  =  (K1*B)  +  (K2*BA2) 

IF  C  <=  0  THEN  K1=K1YNEG:K2=K2YNEG 
:K2=K2YPOS 

Y(J)  =  (K1*C)  +  (K2*CA2) 

IF  D  <=  0  THEN  K1=K1PMNEG:K2=K2PMNEG 


ELSE  Kl=KlLPOS 


ELSE  Kl=KlDPOS 


ELSE  K1=K1YP0S 


PMPOS:K2=K2PMPOS 

PM(J)  =  (K1*D) 
IF  E  <=  0  THEN 

YMPOS:K2=K2YMPOS 

YM(J)  =  (K1*E) 
IF  F  <=  0  THEN 

RMPOS : K2=K2RMPOS 

RM(J)  =  (K1*F) 

NEXT  J 


+  (K2*DA2) 

K1=K1YMNEG : K2=K2 YMNEG 

+  (K2*EA2) 

K1=K1RMNEG: K2=K2RMNEG 

+  (K2*FA2) 


ELSE  K1=K1 


ELSE  K1=K1 


ELSE  K1=K1 
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680   ' PERFORM  INTERACTON  CORRECTIONS 

690   FOR  K  =  1  TO  NOD% 

695   i FIRST  CONSTANT  DETERMINATION 

700     IF  L(K)  >  0  THEN  GOTO  710  ELSE  GOTO  720 

710     DD1=DDDL1P : DD6=DDDL2P : DY1=DYDL1P : DY6=DYDL2  P : DPM1=DP 

MDL1P : DPM6=DPMDL2P : DYM1=DYMDL1P : DYM6=DYMDL2P : DRM1=D 

RMDL1P: DRM6=DRMDL2P 
715     GOTO  73  0 

7  20     DD1=DDDL1N : DD6=DDDL2N : DY1=DYDL1N : DY6=DYDL2N : DPM1=DP 

MDL1N : DPM6=DPMDL2N : DYM1=DYMDL1N : DYM6=DYMDL2N : DRM1=D 

RMDL1N : DRM6=DRMDL2N 
730     IF  D(K)  >  0  THEN  GOTO  740  ELSE  GOTO  750 
740     DL1=DLDD1P : DL6=DLDD2P : DY2=DYDD1P : DY7=DYDD2P : DPM2=DP 

MDD1P : DPM7=DPMDD2  P : DYM2=DYMDD1P : DYM7=DYMDD2  P : DRM2=D 

RMDD1P : DRM7=DRMDD2P 
745     GOTO  760 
750     DL1=DLDD1N : DL6=DLDD2N : DY2=DYDD1N : DY7=DYDD2N : DPM2=DP 

MDD1N : DPM7=DPMDD2N : DYM2=DYMDD1N : DYM7=DYMDD2N : DRM2=D 

RMDD1N : DRM7=DRMDD2N 
760     IF  Y(K)  >  0  THEN  GOTO  770  ELSE  GOTO  780 
770     DL2=DLDY1P: DL7=DLDY2P: DD2=DDDY1P: DD7=DDDY2P: DPM3=DP 

MDY1P : DPM8=DPMDY2P: DYM3=DYMDY1P : DYM8=DYMDY2P : DRM3=D 

RMDY1P: DRM8=DRMDY2P 
775     GOTO  790 
78  0     DL2=DLDY1N : DL7=DLDY2N : DD2=DDDY1N : DD7=DDDY2N : DPM3=DP 

MDY1N : DPM8=DPMDY2N : DYM3=DYMDY1N : DYM8=DYMDY2N : DRM3=D 

RMDY1N: DRM8=DRMDY2N 
790     IF  PM(K)  >  0  THEN  GOTO  800  ELSE  GOTO  810 
800     DL3=DLDPM1P: DL8=DLDPM1P: DD3=DDDPM1P: DD8=DDDPM2P: DY3 

=DYDPM1P : DY8=DYDPM2P : DYM4=DYMDPM1P : DYM9=DYMDPM2P : DR 

M4=DRMDPM1P: DRM9=DRMDPM2P 
805 .    GOTO  820 
810     DL3=DLDPM1N : DL8=DLDPM1N : DD3=DDDPM1N : DD8=DDDPM2N : DY3 

=DYDPM1N : DY8=DYDPM2N : DYM4=DYMDPM1N : DYM9=DYMDPM2N : DR 

M4=DRMDPM1N: DRM9=DRMDPM2N 
820     IF  YM(K)  >  0  THEN  GOTO  830  ELSE  GOTO  840 

8  3  0     DL4=DLDYM1P : DL9=DLDYM2P : DD4=DDDYM1P : DD9=DDDYM2P : DY4 

=DYDYM1P : DY9=DYDYM2P : DPM4=DPMDYM1P : DPM9=DPMDYM2  P : DR 

M5=DRMDYM1P: DRM10=DRMDYM2P 
835     GOTO  850 
84  0     DL4=DLDYM1N: DL9=DLDYM2N: DD4=DDDYM1N: DD9=DDDYM2N: DY4 

=DYDYM1N : DY9=DYDYM2N : DPM4=DPMDYM1N : DPM9=DPMDYM2N : DR 

M5=DRMDYM1N: DRM10=DRMDYM2N 
850     IF  RM(K)  >  0  THEN  GOTO  860  ELSE  GOTO  870 
860     DL5=DLDRM1P: DL10=DLDRM2P: DD5=DDDRM1P: DD10=DDDRM2P: D 

Y5=DYDRM1P : DY10=DYDRM2P : DPM5=DPMDRM1P : DPM10=DPMDRM2 

P : DYM5=DYMDRM1P : DYM10=DYMDRM2  P 
865     GOTO  880 
870     DL5=DLDRM1N:DL10=DLDRM2N:DD5=DDDRM1N:DD10=DDDRM2N:D 

Y5=DYDRM1N:DY10=DYDRM2N:DPM5=DPMDRM1N:DPM10=DPMDRM2 
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N: DYM5=DYMDRM1N : DYM10=DYMDRM2N 

880  » FIRST  INTERACTION  CALCULATION 

890     L1(K)=L(K)-(DL1*D(K) ) - (DL6* (D (K) A2) ) 
900     D1(K)=D(K)-(DD1*L(K) ) - (DD6* (D(K) A2) ) 

950  ' SECOND  CONSTANT  DETERMINATION 

960     IF  L1(K)  >  0  THEN  GOTO  970  ELSE  GOTO  980 

970     DD1=DDDL1P:DD6=DDDL2P:DY1=DYDL1P:DY6=DYDL2P:DPM1=DP 

MDL1P : DPM6=DPMDL2P : DYM1=DYMDL1P : DYM6=DYMDL2  P : DRM1=D 

RMDL1P : DRM6=DRMDL2  P 
975     GOTO  990 
98  0     DD1=DDDL1N : DD6=DDDL2N : DY1=DYDL1N : DY6=DYDL2N : DPM1=DP 

MDL1N : DPM6=DPMDL2N : DYM1=DYMDL1N : DYM6=DYMDL2N : DRM1=D 

RMDL1N : DRM6=DRMDL2N 
990     IF  D1(K)  >  0  THEN  GOTO  1000  ELSE  GOTO  1010 
1000    DL1=DLDD1P:DL6=DLDD2P:DY2=DYDD1P:DY7=DYDD2P:DPM2=DP 

MDD1P : DPM7=DPMDD2P : DYM2=DYMDD1P : DYM7=DYMDD2P : DRM2=D 

RMDD1P: DRM7=DRMDD2P 
1005    GOTO  1020 
1010    DL1=DLDD1N : DL6=DLDD2N : DY2=DYDD1N : DY7=DYDD2N : DPM2=DP 

MDD1N : DPM7=DPMDD2N : DYM2=DYMDD1N : DYM7=DYMDD2N : DRM2=D 

RMDD1N : DRM7=DRMDD2N 
1020    IF  Y1(K)  >  0  THEN  GOTO  1030  ELSE  GOTO  1040 
103  0    DL2=DLDY1P: DL7=DLDY2P: DD2=DDDY1P: DD7=DDDY2P: DPM3=DP 

MDY1P:DPM8=DPMDY2P:DYM3=DYMDY1P:DYM8=DYMDY2P:DRM3=D 

RMDY1P: DRM8=DRMDY2P 
1035    GOTO  1050 
1040    DL2=DLDY1N: DL7=DLDY2N: DD2=DDDY1N: DD7=DDDY2N: DPM3=DP 

MDY1N: DPM8=DPMDY2N: DYM3=DYMDY1N: DYM8=DYMDY2N: DRM3=D 

RMDY1N : DRM8=DRMDY2N 
1050    IF  PM1(K)  >  0  THEN  GOTO  1060  ELSE  GOTO  1080 
1060    DL3=DLDPM1P: DL8=DLDPM1P: DD3=DDDPM1P: DD8=DDDPM2P: DY3 

=DYDPM1P: DY8=DYDPM2P: DYM4=DYMDPM1P: DYM9=DYMDPM2P: DR 

M4=DRMDPM1P: DRM9=DRMDPM2P 
1070    GOTO  1090 
1080    DL3=DLDPM1N : DL8=DLDPM1N : DD3=DDDPM1N : DD8=DDDPM2N : DY3 

=DYDPM1N : DY8=DYDPM2N : DYM4=DYMDPM1N : DYM9=DYMDPM2N : DR 

M4=DRMDPM1N : DRM9=DRMDPM2N 
1090    IF  YM1(K)  >  0  THEN  GOTO  1100  ELSE  GOTO  1110 
1100    DL4=DLDYM1P : DL9=DLDYM2P : DD4=DDDYM1P : DD9=DDDYM2  P : DY4 

=DYDYM1P: DY9=DYDYM2P: DPM4=DPMDYM1P: DPM9=DPMDYM2P: DR 

M5=DRMDYM1P: DRM10=DRMDYM2P 
1105    GOTO  1120 
1110    DL4=DLDYM1N: DL9=DLDYM2N: DD4=DDDYM1N: DD9=DDDYM2N: DY4 

=DYDYM1N : DY9=DYDYM2N : DPM4=DPMDYM1N : DPM9=DPMDYM2N : DR 

M5=DRMDYM1N : DRM10=DRMDYM2N 
1120    IF  RM1(K)  >  0  THEN  GOTO  1130  ELSE  GOTO  1140 
113  0    DL5=DLDRM1P : DL10=DLDRM2P : DD5=DDDRM1P : DD10=DDDRM2P : D 

Y5=DYDRM1P: DY10=DYDRM2P: DPM5=DPMDRM1P: DPM10=DPMDRM2 

P:DYM5=DYMDRM1P:DYM10=DYMDRM2P 
1135    GOTO  1150 

Figure  A. 11  RED.BAS  (cont.) 
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1140    DL5=DLDRM1N : DL10=DLDRM2N : DD5=DDDRM1N : DD10=DDDRM2N : D 

Y5=DYDRM1N : DY10=DYDRM2N : DPM5=DPMDRM1N : DPM10=DPMDRM2 

N: DYM5=DYMDRM1N: DYM10=DYMDRM2N 

1150  ' SECOND  INTERACTION  CALCULATION 

1160    L2(K)=L(K)-(DL1*D1(K) ) - (DL6* (Dl (K) A2) ) 
1170    D2(K)=D(K)-(DD1*L1(K) ) - (DD6* (Dl (K) A2) ) 

1220  '-- — COMPARE  INTERACTION  CALCULATIONS 

1230    DFL=ABS(L1(K)-L2(K) ) : DFD=ABS (Dl (K) -D2 (K) ) 

1250    IF  DFL  <  .0005  AND  DFD  <  .0005  THEN  GOTO  1330 

1260    IF  DFL  >  .0005  THEN  L1(K)=L2(K) 

1270    IF  DFD  >  .0005  THEN  D1(K)=D2(K) 

1320    GOTO  960 

1330    L(K)=L2(K)  :D(K)=D2(K) 

1340    NEXT  K 

23  3  0  » PERFORM  WEIGHT  TARE  CORRECTIONS 

2  3  35  ' — DRAG  AND  LIFT  NEED  TO  BE  CORRECTED  FOR  AXIS  ORIENT 

ATION — 
2340  FOR  K  =  1  TO  NOD% 

2350    ALPHA  =  AOA(K) * (3 . 141593/180)  'CONVERT  TO  RADIANS 
2360    DF(K)=(((-1*D(K) )- (W*SIN (ALPHA) ) ) *COS (ALPHA) )  +  ((( 

-1*L (K) ) - (W*COS (ALPHA) ) ) *SIN (ALPHA) ) 
2370    LF(K)=(( (-1*L(K) ) - (W*COS ( ALPHA) ) ) *COS ( ALPHA) )  -  ((( 

-1*L(K) )-(W*SIN (ALPHA) ) ) *SIN (ALPHA) ) 
2380    NEXT  K 

2390  ' SAVE  REDUCED  DATA  TO  FILE 

2  3  95  CLS: PRINT  "DATA  REDUCTION  COMPLETE,  SAVING  DATA  TO  "; 

FILE$ 
2  4  00  OPEN  FILE$  FOR  OUTPUT  AS  #2 
2410  WRITE  #2,NOD%,CONFIG$,CONF 

2420  WRITE  #2,Q,W 

2421  WRITE  #2,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 

2422  WRITE  #2,CLL, C ALD , C ALY , CALPM , C ALYM , C ALRM , C ALAO A 
2430     FOR  J  =  1  TO  NOD% 

2450     WRITE  #2 , LF ( J) , DF ( J) , Y ( J) , PM( J) , YM( J) ,RM( J) , AOA( J) 
2460     NEXT  J 

2470  ■ RETURN  TO  MAIN  PROGRAM 

248  0  COMMON  PD$,FD$,DX% 

2490  CHAIN  PD$+"MAIN. BAS" , 2110 

2500  ' END  OF  REDUCTION 
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210 

220 

230 
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390 
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CREATE  A  NEW  FILE" 

EDIT  EXISTING  FILE" 

INCREASE  NUMBER  OF  TEST  POINTS" 
IN  AN  EXISTING  FILE" 


i DATA  EDITOR 

GOSUB  180 

COLOR  15,1,7:  KEY  OFF:  CLS 

LOCATE  5,34:  PRINT  "*OPTIONS  MENU*" 

PRINT 

PRINT 

PRINT  TAB (30)  "1. 

PRINT 

PRINT  TAB(30)  "2. 

PRINT 
PRINT  TAB (30)  "3. 
PRINT  TAB (30)  " 
PRINT 

PRINT  TAB (30)  "4.   EXIT  DATA  EDITOR" 

LOCATE  20,5:  INPUT  "ENTER  YOUR  CHOICE  (1, 2 , 3 , 4) " ;REP% 
ON  REP%  GOSUB  390,280,2530,480 
GOTO  3  0 

• INITIALIZING  CONSTANTS ; 

OPTION  BASE  1:  KEY  OFF 

DIM  LF(IOO) ,DF(100) ,YF(100) ,PM(100) ,RM(100) 

DIM  YM(IOO) ,AOA( 100) , COMMAND (10) 

BLANK2$  =  "  ." 

BLANK1$  =  "         " 

BLANK$   =  "  " 

COMMAND$(l)  =  "c":  COMMAND$(2)  =  "r" 

=  "d" 

COMMAND$(4)  =  "e":  COMMANDS (5)  =  "s" 

"u":  COMMAND$ (7)  =  "q" 


COMMANDS (3) 
COMMANDS (6)  = 


RETURN 
i ___ _ _ 


EDIT  A  FILE 

CLS:  COLOR  15,1,7 

LOCATE  5,10:  INPUT  "ENTER  YOUR  FILE  NAME"; FILES 

FILES  =  FD$+FILE$ 

GOSUB  2  680    '  READ  IN  FILE 

PAGE%  =  1:  COL%  =  1:  SAVED%  =  0:  OLDXPOS%  =  0: 

OLDYPOS%  =  0 

GOSUB  1050   'PRINT  BACKGROUND  LINES 

ROW%  =  3:  GOSUB  1740  :   GOSUB  1900      'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 

LOCATE  23,30:  PRINT  "LISTING":  BEEP 

GOSUB  1430   'PRINT  COLUMNS  1-7  AND  Q  VALUE 

LOCATE  23,20:  PRINT  BLANK2$ 

GOTO  470 

i CREATE  A  FILE 

COLOR  15,1,7:  CLS 

LOCATE  5,10:  INPUT  "ENTER  YOUR  FILE  NAME"; FILES 

FILES  =  FD$  +  FILES 

LOCATE  7,10:  INPUT  "ENTER  TEST  CONFIGURATION" ;CONFIG$ 

LOCATE  9,10:  INPUT  "ENTER  THE  NUMBER  OF  TEST  POINTS 

PER  RUN";NOD% 

Figure  A. 12   ADATA.BAS  -  On-Screen  Data  Editor 
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PAGE%  =  1:  COL%  =  1:  SAVED%  =  0:  OLDXPOS%  =  0: 

OLDYPOS%  =  0 

GOSUB  1050   'PRINT  BACKGROUND  LINES 

ROW%  =  3:  GOSUB  1740:  GOSUB  1900       'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 

GOSUB  490   '  INSERT  NEW  DATA 

COMMON  FD$,PD$,DX% 

CHAIN  PD$+"MAIN.BAS",2110 

RETURN 

' ■ CHANGE  OR  INSERT  DATA 

LOCATE  23,12:  PRINT  BLANK$ 

LOCATE  23,12:  INPUT  REP$ :  IF  REP$  =  ""  THEN  GOTO  530 
GOTO  570 

LOCATE  2  3 , 12 :  PRINT  BLANK$ :  LOCATE  23 , 12  :  PRINT 
"INVALID  COMMAND":  BEEP 
FOR  I  =  1  TO  500  STEP  1 
NEXT  I 
GOTO  500 

REP2$  =  RIGHT$(REP$,1) 
Tl%  =  ASC(REP2$) 

IF  Tl%  >=  43  AND  Tl%  <=  57  THEN  GOTO  650 
Tl%  =  Tl%  OR  32:  REP2$  =  CHR$(T1%) 
GOSUB  2  000 

IF  VALID%  =  0  THEN  GOTO  500  ELSE  VALID%  =  0 
ON  FLAG%  GOSUB  1780,1670,2080,2340,2770,2170,2260 
IF  REP2$  <>  "e"  THEN  GOTO  500  ELSE  GOTO  1040 
WHILE  REP2$  <>  "e" 

NEWDATA  =  VAL(REP$) 

IF  QUIT%  =  1  OR  OLDXPOS%  =  0  THEN  QUIT%  =  0:  GOTO 
700    ELSE 

LOCATE  OLDXPOS % , OLDYPOS % , 0 
PRINT  USING  "####. ##";OLDDATA 

LOCATE  23,12:  PRINT  BLANK$ 

R%  =  ROW%  -  2 

IF  COL%  =  1  THEN  LF (R%+18* (PAGE%-1) 
YPOS%  =  7:  GOTO  790:  ELSE 

IF  COL%  =  2  THEN  DF (R%+18* (PAGE%-1) 

=  18:  GOTO  79  0:  ELSE 

IF  COL%  =  3  THEN  YF (R%+18* (PAGE%-1) 

=  29:  GOTO  790:  ELSE 

IF  COL%  =  4  THEN  PM(R%+18* (PAGE%-1) 

=  40:  GOTO  790:  ELSE 

IF  COL%  =  5  THEN  RM (R%+18* (PAGE%-1) 

=  51:  GOTO  790:  ELS 

IF  COL%  =  6  THEN  YM(R%+18* (PAGE%-1) 

=  62:  GOTO  790:  ELS 

IF  COL%  =7  THEN  AOA(R%+18* (PAGE%-1) 

=  71:  GOTO  790:  ELS 

XPOS%  =  ROW% 

LOCATE  XPOS%,YPOS%:  COLOR  0,10 


=  NEWDATA: 


=NEWDATA:  YPOS= 


=NEWDATA:  YPOS: 


=NEWDATA:  YPOS: 


=NEWDATA:  YPOS^ 


=NEWDATA:  YPOS^ 


=NEWDATA:  YPOS: 
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810 
820 
830 

840 
850 
860 
870 
880 
890 
900 
910 
920 
930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 


=  1 


PRINT  USING  ■'####.##»;  NEWDATA:  MODIFIED% 

SAVED%  =  0:  COLOR  15,1,7 

OLDYPOS%  =  YPOS%:  OLDXPOS%  =  XPOS%:  OLDDATA  = 

NEWDATA 

FWD%  =  1 

1%  =  R0W%+1:  GOSUB  1690   'HIGHLIGHT  NEW  ROW 

IF  1%  >  NOD%+2  THEN  FWD%  =  1  ELSE  FWD%  =  0 

IF  1%  >  20  THEN  FWD%  =  1  ELSE  FWD%  =  0 

GOSUB  1800   'HIGHLIGHT  NEW  ROW 

COLOR  15,1,7 

LOCATE  23,12:  PRINT  BLANK$ 

LOCATE  23, 12: INPUT  REP$ :  IF  REP$  = 

GOTO  9  50 

LOCATE  23,12:  PRINT  BLANK$ : 

"INVALID  COMMAND":  BEEP 

GOTO  900 

REP2$  =  RIGHT$(REP$,1) 

Tl%  =  ASC(REP2$) 

IF  Tl%  >=  43  AND  Tl%  <=  57  THEN  GOTO  650 

Tl%  =  Tl%  OR  32:  REP2$  =  CHR$(T1%) 

GOSUB  2  000   'CHECK  COMMANDS 

IF  VALID%  =  0  THEN  GOTO  500  ELSE  VALID%  =  0 

ON  FLAG%  GOSUB  1780,1670,2080,2340,2770,2170,2260 

'COMMAND  SUBROUTINES 

IF  REP2$  <>  "e"  THEN  GOTO  900  ELSE  GOTO  1030 


""  THEN  GOTO  93  0 


LOCATE  23,12:  PRINT 


=  0 


WEND 
RETURN 


2,1 

=  1  TO  8  0 

1,9:  PRINT 


COLOR  15,1,7 

CLS 

LOCATE 

FOR  1% 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

LOCATE 

FOR  1% 

FOR  1% 


■PRINT  BACKGROUND  LINES' 


1,20 

1,32 

1,41 

1,53 

1,62 

1,74 

21,1 

=  1  TO 

=  1  TO 
LOCATE  I 
LOCATE  I 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

80  : 
22 

;,1: 
;,80 


PRINT  CHR$(220) 
"LIFT" 
"DRAG" 
"YAW" 
"PITCH" 
"ROLL" 
"YAW  M." 
"AOA" 


NEXT 


NEXT 

FOR  1%  =  1  TO  21 
LOCATE  I%,14 
LOCATE  I%,25 
LOCATE  I%,36 
LOCATE  I%,47 


PRINT  CHR$(220);:  NEXT 

PRINT  CHR$(22  2) 
PRINT  CHR$(2  22) 


PRINT  CHR$(179) 
PRINT  CHR$(179) 
PRINT  CHR$(179) 
PRINT  CHR$(17  9) 
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1280       LOCATE  I%,58:  PRINT  CHR$(179) 
1290       LOCATE  I%,69:  PRINT  CHR$(179) 

13  00  NEXT 

1310  FOR  1%  =  1  TO  18 

1320      LOCATE  I%+2,2 

1330      PRINT  USING  "###" ;I%+18* (PAGE%-1) ; :  PRINT  » : " 

1340      NEXT 

1350  COLOR  15,4 

1360  LOCATE  22,15:  PRINT  "C  COLUMN   R  ROW   D  DOWN   E  EXIT" 

1370  LOCATE  22,46:  PRINT  "   S  SAVE     U  UP    Q  TUNNEL 

SPEED" 
1380  COLOR  14,0:  LOCATE  2  3,44:  PRINT  CONFIG$ 
1390  COLOR  15,1,7 

14  00  LOCATE  22,2:  PRINT  "[Q  ='» :  LOCATE  22,12:  PRINT  "  ] " 
1410  LOCATE  23,3:  PRINT  "COMMAND  :" 

1420  RETURN 

1430  ' PRINT  COLUMNS  1-7  AND  Q  VALUE 

1440  FOR  1%  =  1  TO  18 

1450     LOCATE  I%+2,7:  PRINT  USING  "####.##" ;LF (I%+18* (PAG 

E%-1) ) 
1460     NEXT 
1470  FOR  1%  =  1  TO  18 
1480     LOCATE  I%+2 , 18 :  PRINT  USING  "####.##" ;DF (I%+18* (PA 

GE%-1) ) 
14  90     NEXT 
1500  FOR  1%  =  1  TO  18 
1510     LOCATE  I%+2,29:  PRINT  USING  "####.##" ;YF (I%+18* (PA 

GE%-1) ) 
1520     NEXT 
1530  FOR  1%  =  1  TO  18 
1540     LOCATE  I%+2,40:  PRINT  USING  "####.##" ;PM(I%+18* (PA 

GE%-1) ) 
1550     NEXT 
1560  FOR  1%  =  1  TO  18 
1570     LOCATE  I%+2,51:  PRINT  USING  "####.##" ;RM(I%+18* (PA 

GE%-1) ) 
1580     NEXT 
1590  FOR  1%  =  1  TO  18 
1600     LOCATE  I%+2,62:  PRINT  USING  "####.##" ;YM(I%+18* (PA 

GE%-1) ) 
1610     NEXT 
1620  FOR  1%  =  1  TO  18 
1630     LOCATE  I%+2,71:  PRINT  USING  "####.##" ;AOA(I%+18* (P 

AGE%-1) ) 
1640     NEXT 

1650  LOCATE  22, 7: PRINT  USING  "###";Q 
16  60  RETURN 
1670  ' ROW  INDEXING 

168  0  FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

169  0  OLDROW%  =  ROW% :  ROW%  ■  ROW%  +  FWD% 
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IF  ROW%  >  NOD%+2  THEN  ROW%  =  3 

IF  ROW%  >  20  THEN  ROW%  =  3 

IF  ROW%  <  3  THEN  ROW%  =  3 

LOCATE  OLDROW%/2:  PRINT  USING  "###" ;OLDROW%-2+18* (PAG 

E%-1) ;:  PRINT  ":" 

COLOR  14,4:  LOCATE  ROW% , 2 

PRINT  USING  "###" ;ROW%-2+18*(PAGE%-l) ; :  PRINT  " : " 

COLOR  15,1,7 

RETURN 

i COLUMN  INDEXING 

FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

OLDCOL%  =  COL%:  COL%  =  COL%  +  FWD% 

IF  COL%  >  7  THEN  COL%  =  1 

IF  COL%  <  1  THEN  COL%  =  7 

IF  OLDCOL%  =  1  THEN  LOCATE  1,9:  PRINT  "LIFT":  GOTO 

1900  ELSE 

IF  OLDCOL%  =  2  THEN  LOCATE  1,20:  PRINT 

1900  ELSE 

IF  OLDCOL%  =  3  THEN  LOCATE  1,32:  PRINT 

1900  ELSE 

IF  OLDCOL%  =  4  THEN  LOCATE  1,41:  PRINT 

1900  ELSE 

IF  OLDCOL%  =  5  THEN  LOCATE  1,53:  PRINT 

1900  ELSE 

IF  OLDCOL%  =  6  THEN  LOCATE  1,62:  PRINT 

1900  ELSE 

IF  OLDCOL%  =  7  THEN  LOCATE  1,74:  PRINT 

1900  ELSE 

COLOR  14,4 

IF  COL%  =  1  THEN  LOCATE  1,9:  PRINT  "LIFT":  GOTO  1980 

ELSE 

IF  COL%  =  2  THEN  LOCATE  1,20 

ELSE 

IF  COL%  =  3  THEN  LOCATE  1,3  2 

ELSE 

IF  COL%  =  4  THEN  LOCATE  1,41 

ELSE 

IF  COL%  =  5  THEN  LOCATE  1,53 

ELSE 

IF  COL%  =  6  THEN  LOCATE  1,62 

1980  ELSE 

IF  COL%  =  7  THEN  LOCATE  1,74 

ELSE 

COLOR  15,1,7 

RETURN 

i CHECK  COMMANDS 

1%  =  0 

WHILE  1%  <  7  AND  VALID%  =  0 

1%  =  1%  +  1 

IF  REP2$  =  COMMAND$(I%)  THEN  FLAG%  =  1%:  VALID%  =  1 


"DRAG":  GOTO 
"YAW" :  GOTO 
"PITCH":  GOTO 
"ROLL":  GOTO 
"YAW  M.":  GOTO 
"AOA":  GOTO 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


"DRAG":  GOTO  198  0 
"YAW":  GOTO  198  0 
"PITCH":  GOTO  1980 
"ROLL":  GOTO  19  80 
"YAW  M.":  GOTO 
"AOA":  GOTO  19  8  0 
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2050     WEND 

2060     IF  VALID%  =  0  THEN  GOSUB  2420   'ERROR  MESSAGE 

2070  RETURN 

2080  ' SCROLL  DOWN 

2  090  FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

2100  PAGE%  =  PAGE%  +  FWD% 

2110  IF  PAGE%  <  1  THEN  PAGE%  =  1 

2120  GOSUB  1050   'PRINT  BACKGROUND  LINES 

213  0  ROW%  =  3:  COL%  =  1:  GOSUB  174  0:  GOSUB  1900  'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 
2140  GOSUB  143  0  'PRINT  COLUMNS  1-7  AND  Q  VALUE 
2150  OLDXPOS%  =  0 
2160  RETURN 

2170  ' SCROLL  UP ■ 

2180  FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

2190  PAGE%  =  PAGE%  -  FWD% 

2200  IF  PAGE%  <  1  THEN  PAGE%  =  1 

2210  GOSUB  1050  'PRINT  BACKGROUND  LINES 

2  220  ROW%  =  3:  COL%  =  1:  GOSUB  174  0:  GOSUB  1900  'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 
22  3  0  GOSUB  143  0  '  PRINT  COLUMNS  1-7  AND  Q  VALUE 
2240  OLDXPOS%  =  0 
2250  RETURN 

2260  ' PRINT  OUT  Q  VALUE 

2270  LOCATE  23,12:  PRINT  BLANK$ 

2280  LOCATE  23,12:  INPUT  "Q  =";Q 

2290  LOCATE  2  3,12:  PRINT  BLANK$ 

2300  COLOR  0,10 

2  310  LOCATE  22,7:  PRINT  USING  »###», -Q 

2320  COLOR  15,1,7 

2330  RETURN 

2340  '- ■ EXIT  EDITOR 

2350  IF  SAVED%  =  1  OR  MODIFIED%  =  0  THEN  GOTO  2410 

2  3  60  LOCATE  23,20:  PRINT  BLANK2$:  BEEP 

2370  LOCATE  23,20:  INPUT  "SAVE  FILE  (Y/N)";REP$ 

2380  IF  REP$  =  "N"  OR  REP$  =  "n"  THEN  GOTO  2410 

2  3  90  GOSUB  2770    'SAVE  FILE 

2400  LOCATE  23,20:  PRINT  BLANK2$ 

2410  RETURN 

2420  ' ■ ERROR  MESSAGE ■ 

243  0  LOCATE  23,20:  PRINT  BLANK2$ 

244  0  LOCATE  23,20:  PRINT  "  INVALID  COMMAND" 
2450  FOR  F  =  300  TO  500  STEP  100 

2460     SOUND  F,2 

2470     SOUND  32767,2 

2480     NEXT 

2490  FOR  I  =  1  TO  500  STEP  1 

2500  NEXT  I 

2510  LOCATE  23,20:  PRINT  BLANK2$ 

252  0  RETURN 
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2530  ' ADD  TEST  POINTS  TO  FILE 

2540  CLS 

2550  LOCATE  3,5:  PRINT"THIS  OPTION  ALLOWS  YOU  TO  INCREASE 

THE  NUMBER  OF  TEST" 
2560  LOCATE  4,5:  PRINT"POINTS  OF  AN  EXISTING  FILE" 
2570  LOCATE  7,5:  INPUT"WHAT  IS  THE  NAME  OF  YOUR  FILE";FILE$ 
2580  FILE$  =  FD$  +  FILE$ 
2590  GOSUB  2680 
2610  LOCATE  10,5:  INPUT"WHAT  IS  THE  NEW  NUMBER  OF  TUNNEL 

RUNS";NOQ% 
2620  LOCATE  13,5:  INPUT"WHAT  IS  THE  NEW  NUMBER  OF  TEST 

POINTS  PER  RUN";NOD% 
2630  GOSUB  2770 
2660  LOCATE  20,5:  PRINT"NOW  YOU  CAN  CALL  UP  YOUR  FILE  AND 

ADD  IN  THE  NEW  POINTS" 
2  67  0  RETURN 

2680  ' READ  IN  FILE 

2  690  OPEN  FILE$  FOR  INPUT  AS  #1 
27  00  INPUT  #l,NOD%,CONFIG$,CONF 
2710  INPUT  #1,Q,W 

2714  INPUT  #l,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 

2715  INPUT  #l,CLL,CALD,CALY,CALPM,CALYM,CALRM,CALAOA 
2720    FOR  J  =  1  TO  NOD% 

2730    INPUT  #1,LF(J) ,DF(J) ,YF(J) ,PM(J) ,RM(J) ,YM(J) ,AOA(J) 
2740    NEXT  J 
2750  CLOSE  #1 
2760  RETURN 

2770  ■ SAVE  FILE 

2780  LOCATE  23,20:  PRINT  BLANK2$:  LOCATE  23,30:  PRINT 

"SAVING  FILE" 
2790  OPEN  FILE$  FOR  OUTPUT  AS  #2 
2  8  00  WRITE  #2,NOD%,CONFIG$,CONF 
2810  WRITE  #2,Q,W 

2814  WRITE  #2,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 
2  815  WRITE  #2,CLL,CALD,CALY,CALPM,CALYM,CALRM,CALAOA 
2820   FOR  J  =  1  TO  NOD% 

2830     WRITE  #2 , LF ( J) , DF (J) , YF (J) , PM( J) , RM( J) , YM( J) , AOA( J) 
2840     NEXT  J 
2850  CLOSE  #2 
2  8  60  RETURN 
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10  REM COMP.BAS;   CALCULATE  Cd,  CI  AND  E.  F.  A.— 

15  REM 

16  DIM  LF(100) ,DF(100) ,YF(100) ,PM(100) ,YM(100) ,RM(100) , 
AOA(IOO) 

17  DIM  CD(IOO) ,CL(100) ,EFA(100) 

2  0  COLOR  15, 1,4: KEY  OFF:  CLS 

3  0  PRINT  "WHAT  IS  THE  NAME  OF  YOUR  FILE  CONTAINING  THE"; 
35  PRINT  "  CALIBRATED" 

40  INPUT  "WIND  TUNNEL  DATA";FILE$ 

50  FILE$  =  FD$+  FILE$ 

60  GOSUB  680         'READ  IN  DATA  FROM  FILE 

7  0  REM  —  FOR  THE  FOLLOWING,  AREAL  IS  MEASURED  IN  SQUARE 

FEET 

90  AREAL  =60/144 

270  FOR  J  =  1  TO  NOD% 
280        CL(J)  =  LF(J)/(Q*AREAL) 
290        CD (J)  =  DF(J)/(Q*AREAL) 
295        EFA(J)  =LF(J)/Q 
3  00        NEXT  J 

310  ' STORE  DATA  FOR  Cd  VS.  CI 

311  INPUT"WHAT  IS  THE  NAME  FOR  THE  FILE  TO  STORE  Cd  VS. 
CI  DATA";N$ 

312  Y$="Cd":X$="Cl" 

32  0  OPEN  FD$+N$  FOR  OUTPUT  AS  #1. 

330  WRITE  #l,NOD%,Y$,X$,CONFIG$,CONF 

340  WRITE  #1,Q 

3  50        FOR  J  =   1  TO  NOD% 

360  WRITE  #1,CD(J) ,CL(J) 

370  NEXT  J 

380  CLOSE  #1 

390  • STORE  DATA  FOR  Cd  VS.  C1*C1 

391  INPUT"WHAT  IS  THE  NAME  FOR  THE  FILE  TO  STORE  Cd  VS. 
C1A2  DATA";N$ 

392  Y$="Cd":X$="ClA2" 

4  00  OPEN  FD$+N$  FOR  OUTPUT  AS  #2 

410  WRITE  #2,NOD%,Y$,X$,CONFIG$,CONF 

42  0  WRITE  #2,Q 

43  0        FOR  J  =   1  TO  NOD% 
440  B  =  CL(J)*CL(J) 
450           WRITE  #2,CD(J) ,B 
4  60           NEXT  J 

470  CLOSE  #2 

480  ' STORE  DATA  FOR  Cl  VS.  AOA 

481  INPUT"WHAT  IS  THE  NAME  FOR  THE  FILE  TO  STORE  THE  Cl 
VS.  AOA  DATA";N$ 

482  Y$="C1":X$="A0A" 

490  OPEN  FD$+N$  FOR  OUTPUT  AS  #2 

500  WRITE  #2,NOD%,Y$,X$,CONFIG$,CONF 

510  WRITE  #2,Q 

52  0        FOR  J  =   1  TO  NOD% 

Figure  A. 13   COMP.BAS  -  Program  to  Calculate  Parameters 
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530  WRITE  #2,CL(J) ,AOA(J) 

540  NEXT  J 

550  CLOSE  #2 

560  ' STORE  DATA  FOR  E.F.A.  VS  AOA 

561  INPUT "WHAT  IS  THE  NAME  FOR  THE  FILE  TO  STORE  THE 

.F.A.  VS.  AOA  DATA";N$ 

562  Y$="E.F.A.":X$="AOA" 

570  OPEN  FD$+N$  FOR  OUTPUT  AS  #1 

580  WRITE  #l,NOD%,Y$,X$,CONFIG$,CONF 

590  WRITE  #1,Q 

600        FOR  J  =   1  TO  NOD% 

610  WRITE  #1,EFA(J) ,AOA(J) 

62  0  NEXT  J 

630  CLOSE  #1 

640  ' RETURN  TO  MAIN  PROGRAM 

650  COMMON  FD$,PD$,DX% 

660  CHAIN  PD$+"MAIN.BAS",2110 

670  • READ  IN  DATA  FROM  CALIBRATED  FILE 

68  0  OPEN  FILE$  FOR  INPUT  AS  #2 

69  0  INPUT  #2,NOD%,CONFIG$,CONF 

700  INPUT  #2,Q,W 

701  INPUT  #2,ZL,ZD,ZY,ZPM,ZYM,ZRM,ZAOA 

7  02  INPUT  #  2 , CLL , CALD , CALY , CALPM , CALYM , CALRM , CALAOA 

710        FOR  J  =  1  TO  NOD% 

720        INPUT  #2,LF(J) ,DF(J) ,YF(J) ,PM(J) ,RM(J) ,YM(J) , 

AOA (J) 
730        NEXT  J 

740  CLOSE  #2 

7  50  RETURN 
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CREATE  A  NEW  FILE" 

EDIT  EXISTING  FILE" 

INCREASE  NUMBER  OF  TEST  POINTS" 
IN  AN  EXISTING  FILE" 


H 

U3 


' ■ -DATA  EDITOR ■— 

GOSUB  180 

COLOR  15,1,7:  KEY  OFF:  CLS 

LOCATE  5,34:  PRINT  "*OPTIONS  MENU*" 

PRINT 

PRINT 

PRINT  TAB (30)  "1. 

PRINT 

PRINT  TAB (30)  "2. 

PRINT 
PRINT  TAB (3  0)  "3. 
PRINT  TAB (30)  " 
PRINT 

PRINT  TAB (30)  "4.   EXIT  DATA  EDITOR" 

LOCATE  20,5:  INPUT  "ENTER  YOUR  CHOICE  (1 , 2 , 3 , 4) " ;REP% 
ON  REP%  GOSUB  380,270,2200,490 
GOTO  3  0 

i INITIALIZING  CONSTANTS 

OPTION  BASE  1:  KEY  OFF 

DIM  Y(100) ,X(100) ,COMMAND(10) 

BLANK2$  =  " 

BLANK1$  =  "         " 

BLANK$   =  "  " 

COMMAND$(l)  =  "c":COMMAND$(2)  =  "r" :  COMMAND$(3)  =  "d" 

COMMAND$(4)  =  "e":  COMMAND$(5)  =  "s":  COMMAND$(6)  = 

"u":  COMMAND$(7)  =  "q" 

RETURN 

i  — EDIT  A  FILE 

CLS:  COLOR  15,1,7 

LOCATE  5,10:  INPUT  "ENTER  YOUR  FILE  NAME";FILE$ 

FILE$  =  FD$  +  FILE$ 

GOSUB  2  3  50    ■  READ  IN  FILE 

PAGE%  =  1:  COL%  =  1:  SAVED%  =  0:  OLDXPOS%  =  0: 

OLDYPOS%  =  0 

GOSUB  1010   'PRINT  BACKGROUND  LINES 

ROW%  =  3:  GOSUB  152  0  :   GOSUB  162  0      'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 

LOCATE  23,30:  PRINT  "LISTING":  BEEP 

GOSUB  1360   'PRINT  COLUMNS  1-7  AND  Q  VALUE 

LOCATE  23,20:  PRINT  BLANK2$ 

GOTO  48  0 

' CREATE  A  FILE 

COLOR  15,1,7:  CLS 

LOCATE  5,10:  INPUT  "ENTER  YOUR  FILE  NAME";FILE$ 

FILR$  =  FD$  +  FILE$ 

LOCATE  7,10:  INPUT  "ENTER  TEST  CONFIGURATION" ;CONFIG$ 

LOCATE  9,10:  INPUT  "ENTER  THE  NUMBER  OF  TEST  POINTS 

PER  RUN";NOD% 

LOCATE  11,10:  INPUT  "ENTER  TITLE  FOR  X-VALUES  (i.e. 

CL,  A0A)";X1$ 

Figure  A. 14   BDATA.BAS   -  On-Screen  Data  Editor 
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440   LOCATE  13,10:  INPUT  "ENTER  TITLE  FOR  Y-VALUES  (i.e. 

Cd,  C1)";Y1$ 
450   PAGE%  =  1:  COL%  =  1:  SAVED%  =  0:  OLDXPOS%  =  0: 

OLDYPOS%  =  0 
4  60   GOSUB  1010   'PRINT  BACKGROUND  LINES 
47  0   ROW%  =  3:  GOSUB  1520:  GOSUB  1620       'HIGHLIGHT  FIRST 

ROW  AND  COLUMN 
480   GOSUB  500   '  INSERT  NEW  DATA 

490  COMMON  FD$,  PD$ , DX% 

491  CHAIN  PD$+"MAIN.BAS",2110 

492  RETURN 

500   ■ CHANGE  OR  INSERT  DATA 

510   LOCATE  23,12:  PRINT  BLANK$ 

52  0   LOCATE  23,12:  INPUT  REP$ :  IF  REP$  =  ""  THEN  GOTO  540 

530   GOTO  580 

540   LOCATE  23,12:  PRINT  BLANK$ :  LOCATE  23,12:  PRINT 

"INVALID  COMMAND":  BEEP 
550    FOR  I  =  1  TO  500  STEP  1 
560    NEXT  I 
570   GOTO  510 

580   REP2$  =  RIGHT$(REP$,1) 
590   Tl%  =  ASC(REP2$) 

600   IF  Tl%  >=  43  AND  Tl%  <=  57  THEN  GOTO  660 
610   Tl%  =  Tl%  OR  32:  REP2$  =  CHR$(T1%) 
620   GOSUB  1670 

630   IF  VALID%  =  0  THEN  GOTO  510  ELSE  VALID%  =  0 
640   ON  FLAG%  GOSUB  1550,1450,1750,2010,2440,1840,1930 
650   IF  REP2$  <>  "e"  THEN  GOTO  510  ELSE  GOTO  1000 
660   WHILE  REP2$  <>  "e" 
670      NEWDATA  =  VAL(REP$) 
680      IF  QUIT%  =  1  OR  OLDXPOS%  =  0  THEN  QUIT%  =  0:  GOTO 

710    ELSE 
690  LOCATE  OLDXPOS% ,OLDYPOS% , 0 

7  00  PRINT  USING  "####.##" ;OLDDATA 

710      LOCATE  23,12:  PRINT  BLANK$ 
72  0      R%  =  ROW%  -  2 
730      IF  COL%  =  1  THEN  Y (R%+18* (PAGE%-1) ) =NEWDATA:  YPOS% 

=  18:  GOTO  750:  ELS 
740      IF  COL%  =  2  THEN  X(R%+18* (PAGE%-1) ) =NEWDATA:  YPOS% 

=  43:  GOTO  750:  ELSE 
750      XPOS%  =  ROW% 

7  60      LOCATE  XPOS%,YPOS%:  COLOR  0,10 

770      PRINT  USING  »####.##»;  NEWDATA:  MODIFIED%  =  1 
78  0      SAVED%  =  0:  COLOR  15,1 
790      OLDYPOS%  =  YPOS%:  OLDXPOS%  =  XPOS%:  OLDDATA  = 

NEWDATA 
800   FWD%  =  1 
810   1%  =  R0W%+1:  GOSUB  147  0   'HIGHLIGHT  NEW  ROW 

82  0   IF  1%  >  NOD%+2  THEN  FWD%  =  1  ELSE  FWD%  =  0 

83  0   IF  1%  >  20  THEN  FWD%  =  1  ELSE  FWD%  =  0 
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GOSUB  1570   'HIGHLIGHT  NEW  ROW 
COLOR  15,1,7 

LOCATE  23,12:  PRINT  BLANK$ 

LOCATE  23,12:  INPUT  REP$ :  IF 

890 

GOTO  910 

LOCATE  23,12 


REP$  =  »"  THEN  GOTO 


PRINT  BLANK$: 
BEEP 


LOCATE  2  3,12:  PRINT 
"INVALID  COMMAND": 
GOTO  860 

REP2$  =  RIGHT$(REP$,1) 
Tl%  =  ASC(REP2$) 

IF  Tl%  >=  43  AND  Tl%  <=  57  THEN  GOTO  660 
Tl%  =  Tl%  OR  32:  REP2$  =  CHR$(T1%) 
GOSUB  i670   'CHECK  COMMANDS 

IF  VALID%  =  0  THEN  GOTO  510  ELSE  VALID%  =  0 
ON  FLAG%  GOSUB  1550,1450,1750,2010,2440,1840,1930 
'COMMAND  SUBROUTINES 

IF  REP2$  <>  "e"  THEN  GOTO  860  ELSE  GOTO  990 
WEND 


RETURN 


■PRINT  BACKGROUND  LINES- 


80  :  PRINT  CHR$(220) 
PRINT  "Y  =  ";Y1$ 
PRINT  "X  =  ";X1$ 
PRINT  "COMMANDS" 

PRINT  CHR$(220) 


NEXT 


80 
22 

I%,1: 
I%,80 


NEXT 


COLOR  15,1,7 
CLS 
LOCATE  2,1 
FOR  1%  =  1  TO 
LOCATE  1,18: 
LOCATE  1,43: 
LOCATE  1,68: 
LOCATE  21,1 
FOR  1%  =  1  TO 
FOR  1%  =  1  TO 
LOCATE 
LOCATE 
NEXT 

FOR  1%  =  1  TO  21 
LOCATE  I%,28 
LOCATE  I%,53 
NEXT 

FOR  1%  =  1  TO  18 
LOCATE  I%+2,2 

PRINT  USING  "###";I%+18*(PAGE%-1) ; :  PRINT 
NEXT 
COLOR  15,4 


PRINT  CHR$(222) 
:  PRINT  CHR$(222) 


PRINT  CHR$(179) 
PRINT  CHR$(179) 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 
LOCATE 


3,66 
4,66 
5,66 
6,66 
7,66 
8,66 
9,66 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
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"Q 
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DOWN 
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SPEED" 
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COLOR  14,0:  LOCATE  23,44:  PRINT  CONFIG$ 

COLOR  15,1,7 

LOCATE  22,2:  PRINT  " [Q  =" :  LOCATE  22,12:  PRINT  "]" 

LOCATE  2  3,3:  PRINT  "COMMAND  :" 

RETURN 

i PRINT  COLUMNS  1-7  AND  Q  VALUE 

FOR  1%  =  1  TO  18 

LOCATE  I%+2,18:  PRINT  USING  "####.##" ;Y (I%+18* (PAG 

E%-1) ) 

NEXT 
FOR  1%  =  1  TO  18 

LOCATE  I%+2,43:  PRINT  USING  "#### . ##" ;X (I%+18* (PAG 

E%-1) ) 

NEXT 
LOCATE  2 2, 7: PRINT  USING  »###";Q 
RETURN 

i ROw  INDEXING 

FWD%  =  VAL(REP$)  :  IF  FWD%  =  0  THEN  FWD%  =•  1 

OLDROW%  =  ROW%:  ROW%  =  ROW%  +  FWD% 

IF  ROW%  >  NOD%+2  THEN  ROW%  =  3 

IF  ROW%  >  20  THEN  ROW%  =  3 

IF  ROW%  <  3  THEN  ROW%  =  3 

LOCATE  OLDROW%,2:  PRINT  USING  "###" ;OLDROW%-2+18* (PAG 

E%-1) ; :  PRINT  ":" 

COLOR  14,4:  LOCATE  ROW% , 2 :  PRINT  USING  "###" ;ROW%-2+l 

8*(PAGE%-1) ; : PRINT  ":" 

COLOR  15,1,7 

RETURN 

i COLUMN  INDEXING 

FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

OLDCOL%  =  COL%:  COL%  =  COL%  +  FWD% 

IF  COL%  >  2  THEN  COL%  =  1 

IF  COL%  <  1  THEN  COL%  =  2 

IF  OLDCOL%  =  1  THEN  LOCATE  1,18:  PRINT  "Y  =  ";Y1$: 

GOTO  1620  ELSE 

IF  OLDCOL%  =  2  THEN  LOCATE  1,43:  PRINT  "X  =  ";X1$: 

GOTO  1620  ELSE 

COLOR  14,4 

IF  COL%  =  1  THEN  LOCATE  1,18:  PRINT  "Y  =  ";Y1$:  GOTO 

1650  ELSE 

IF  COL%  =  2  THEN  LOCATE  1,43:  PRINT  "X  =  ";X1$:  GOTO 

1650  ELSE 

COLOR  15,1,7 

RETURN 

i CHECK  COMMANDS 

1%  =  0 

WHILE  1%  <  7  AND  VALID%  =  0 
1%  =  1%  +  1 

IF  REP2$  =  COMMANDS (1%)  THEN  FLAG%  =  1%:  VALID%  =  1 
WEND 
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173  0     IF  VALID%  =  0  THEN  GOSUB  2090   'ERROR  MESSAGE 
1740  RETURN 

1750  ' SCROLL  DOWN ■ 

1760  FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

1770  PAGE%  =  PAGE%  +  FWD% 

1780  IF  PAGE%  <  1  THEN  PAGE%  =  1 

1790  GOSUB  1010   'PRINT  BACKGROUND  LINES 

1800  ROW%  =  3:  COL%  =  1:  GOSUB  1520:  GOSUB  162  0  'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 
1810  GOSUB  13  60  'PRINT  COLUMNS  1-7  AND  Q  VALUE 

182  0  OLDXPOS%  =  0 

183  0  RETURN 

1840  • SCROLL  UP 

1850  FWD%  =  VAL(REP$) :  IF  FWD%  =  0  THEN  FWD%  =  1 

1860  PAGE%  =  PAGE%  -  FWD% 

1870  IF  PAGE%  <  1  THEN  PAGE%  =  1 

1880  GOSUB  1010  'PRINT  BACKGROUND  LINES 

1890  ROW%  =  3:  COL%  =  1:  GOSUB  1520:  GOSUB  1620  'HIGHLIGHT 

FIRST  ROW  AND  COLUMN 
1900  GOSUB  13  60  '  PRINT  COLUMNS  1-7  AND  Q  VALUE 
1910  OLDXPOS%  =  0 
192  0  RETURN 

1930  ' PRINT  OUT  Q  VALUE >- 

1940  LOCATE  23,12:  PRINT  BLANK$ 

1950  LOCATE  23,12:  INPUT  "Q  =";Q 

1960  LOCATE  23,12:  PRINT  BLANK$ 

1970  COLOR  0,10 

198  0  LOCATE  22,7:  PRINT  USING  "###"' ;Q 

1990  COLOR  15,1,7 

2  000  RETURN 

2010  ' ■ EXIT  EDITOR ■ 

2  02  0  IF  SAVED%  =  1  OR  MODIFIED%  =  0  THEN  GOTO  2  08  0 

2  03  0  LOCATE  23,20:  PRINT  BLANK2$:  BEEP 

2  04  0  LOCATE  23,20:  INPUT  "SAVE  FILE  (Y/N)";REP$ 

2  050  IF  REP$  =  "N"  OR  REP$  =  "n"  THEN  GOTO  2  08  0 

2  060  GOSUB  2440    'SAVE  FILE 

2  070  LOCATE  23,20:  PRINT  BLANK2$ 

2  080  RETURN 

2090  ' ERROR  MESSAGE 

2100  LOCATE  23,20:  PRINT  BLANK2$ 

2110  LOCATE  23,20:  PRINT  "  INVALID  COMMAND" 

2120  FOR  F  =  300  TO  500  STEP  100 

213  0     SOUND  F,2 

2140     SOUND  32767,2 

2150     NEXT 

2160  FOR  I  =  1  TO  500  STEP  1 

217  0  NEXT  I 

218  0  LOCATE  2  3,20:  PRINT  BLANK2$ 
2190  RETURN 

2200  ' ADD  TEST  POINTS  TO  FILE 
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2210  CLS 

2220  LOCATE  3,5:  PRINT"THIS  OPTION  ALLOWS  YOU  TO  INCREASE 

THE  NUMBER  OF  TEST" 
2230  LOCATE  4,5:  PRINT"POINTS  OF  AN  EXISTING  FILE" 
2240  LOCATE  7,5:  INPUT"WHAT  IS  THE  NAME  OF  YOUR  FILE";FILE$ 
2250  FILE$  =  FD$  +  FILE$ 
2260  GOSUB  2350 
2290  LOCATE  13,5:  INPUT"WHAT  IS  THE  NEW  NUMBER  OF  TEST 

POINTS  PER  RUN";NOD% 
2300  GOSUB  2440 
2330  LOCATE  20,5:  PRINT"NOW  YOU  CAN  CALL  UP  YOUR  FILE  AND 

ADD  IN  THE  NEW  POINTS" 
2340  RETURN 

2350  ■ READ  IN  FILE 

23  60  OPEN  FILE$  FOR  INPUT  AS  #1 

2  370  INPUT  #l,NOD%,Yl$,Xl$,CONFIG$,CONF 

2380  INPUT  #1,Q 

2390     FOR  J  =  1  TO  NOD% 

2400        INPUT  #1,Y(J),X(J) 

2410        NEXT  J 

2420  CLOSE  #1 

2430  RETURN 

2440  • SAVE  FILE 

2  450  LOCATE  23,20:  PRINT  BLANK2$:  LOCATE  23,30:  PRINT 

"SAVING  FILE" 
2  460  OPEN  FILE$  FOR  OUTPUT  AS  #2 
2470  WRITE  #2 ,NOD% , Yl$ ,X1$ , CONFIG$ , CONF 
2480  WRITE  #2,Q 
2490     FOR  J  =  1  TO  NOD% 
2500        WRITE  #2,Y(J),X(J) 
2  510        NEXT  J 
2520  CLOSE  #2 
253  0  RETURN 
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PLOTTER- 


XMAX=XMIN:  YMIN=Y(1) 


i , PLOTTING  ROUTINE — ■ — — 

DIM  X(300) ,Y(300) ,8(300) ,C(300) ,D(300) 

KEY  OFF:  SCREEN  2 

HDPOS%  =  44:  B$  =  "  " 

GOSUB  2  60 

OPT$  =  CURVE$(1) 

ERASE  X,Y,FILESTK$,CURVE$,B,C,D 

CHAIN  PD$+"MAIN.BAS",2110,ALL 

i- read  IN  FILE ■- 

OPEN  FILE$  FOR  INPUT  AS  #2 
INPUT  #2,N0D%,Y1$,X1$,H$,C0NF 
INPUT  #2,Q 
FOR  1%  =  1  TO  NOD% 

INPUT  #2,Y(I%) ,X(I%) 

NEXT 
CLOSE  #2 
RETURN 

i INITIALIZING 

FOR  L%  =  1  TO  NOF% 

FILE$  =  FILESTK$(L%) 

GOSUB  160 

IF  L%=1  THEN  XMIN=X(1)- : 
YMAX=YMIN 

FOR  I%=1  TO  NOD% 

IF  X(I%)<XMIN  THEN  XMIN=X(I%) 
XMAX  THEN  XMAX=X(I%) 

IF  Y(I%)<YMIN  THEN  YMIN=Y(I%) 
YMAX  THEN  YMAX=Y(I%) 
NEXT  1% 

NEXT  L% 
CLS 

PRINT  "AFTER  CURVES  ARE  PLOTTED  YOU  HAVE  THE  OPTION 
OF  OBTAINING  A  HARD" 

PRINT  "COPY  BY  PRESSING  CTRL-PRTSC.   THE  CURVE  WILL 
REMAIN  ON  THE  SCREEN" 
PRINT  "UNTIL  YOU  PRESS  THE  F2  KEY." 
PRINT 

INPUT  "HEADDING  FOR  PLOT  =";HEAD$ 

T1%=LEN(HEAD$) :  IF  Tl%>60  THEN  HEAD$=LEFT$ (HEAD$ , 60) 
INPUT  "X  TITLE  FOR  THE  PLOT  =";XTITLE$ 

T1%=LEN(XTITLE$) :  IF  Tl%>60  THEN  XTITLE$=LEFT$ (XTITLE 
$,60) 

INPUT  "Y  TITLE  FOR  THE  PLOT  =";YTITLE$ 

T1%=LEN(YTITLE$) :  IF  Tl%>60  THEN  YTITLE$=LEFT$ (YTITLE 
$,60) 

PRINT  "MINIMUM  X  VALUE  =  [ DEFAULT : " ; XMIN ; " ] " ; :  INPUT 
Tl$ 

IF  Tl$  <>  ""  THEN  TEMP=VAL(T1$)  ELSE  TEMP=XMIN 
IF  TEMP<=XMIN  THEN  0K%=1  ELSE  OK%=0 
IF  Tl$  <>  ""  THEN  TEMP3  =XMIN:  XMIN  =  TEMP 


ELSE  IF 
ELSE  IF 


X(I%)> 
Y(I%)> 
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460  PRINT  "MAXIMUM  X  VALUE  =  [DEFAULT: " ;XMAX; "]"; :  INPUT 

Tl$ 

470  IF  Tl$  <>  ""  THEN  XMAX=VAL ( T 1 $ ) 

480  PRINT  "MINIMUM  Y  VALUE  =  [ DEFAULT :"; YMIN ;"]"; :  INPUT 

Tl$ 

490  IF  Tl$  <>  ""  THEN  YMIN=VAL(T1$) 

500  PRINT  "MAXIMUM  Y  VALUE  =  [DEFAULT: " ;YMAX; "]"; :  INPUT 

Tl$ 

510  IF  Tl$  <>  ""  THEN  YMAX=VAL(T1$) 

52  0  RANGE=XMAX-XMIN 

53  0  IF  0K%=1  THEN  PTR%=1:  STKPTR%=1:  PGSTK%(1)=1  ELSE 

GOSUB  580 

54  0  CLS:  SCREEN  2 
550  GOSUB  840 
560  RETURN 

570  ' SET  PAGES 

580  TEMP2=(TEMP-TEMP3)/RANGE+1:  PG%=INT (TEMP2) 

59  0  STKPTR%=PG%+1 

600  FOR  1%  =  1  TO  STKPTR% 

610       T2=TEMP-RANGE*(I%-1) 

620       J%=1 

630       WHILE  X(J%)<T2 

640  J%=J%+1 

650  WEND 

660       PGSTK%(PG%-I%+2)=J% 

670  NEXT  1% 

680  PTR%=PGSTK%(STKPTR%) 

690  RETURN 

700  • PLOT  WITHOUT  GRID 

710  LINE (ORIGINX% , ORIGINY%) - (ORIGINX%+7*INTERVALX% , TEMP4% 
)fl*B 

72  0  X%=ORIGINX% 

73  0  FOR  1%  =  1  TO  8 

740       LINE  (X%,162)-(X%,165) 

750       X%=X%+INTERVALX% 

760       NEXT 

770  Tl%=ORIGINX%-8:  T2%=ORIGINY% 

780  FOR  1%  =  1  TO  5 

790       LINE  (Tl%,T2%)-(ORIGINX%,T2%) 

8  00       T2%=T2%+INTERVALY% 

810       NEXT 

82  0  RETURN 

330  i PLOTTING  ROUTINES 

840  ORIGINY%=12:  ORIGINX%=72:  INTERVALY%=3  0 :  INTERVALX%=80 

850  TEMP4%=ORIGINY%+5*INTERVALY% 

870  VIEW(0,0)-(639,163) :  CLS:  VIEW 

890  GOSUB  710 

900  IF  LEFT$(GRD$,1)="G"  THEN  GOSUB  2400 

910  X%=ORIGINX%:  Y%=ORIGINY% 

92  0  FOR  L%=1  TO  NOF% 

Figure  A. 15   PLOTTER. BAS  (cont.) 
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IF  N0F%=1  THEN  GOTO  960 

FILE$=FILESTK$ (L%) 

GOSUB  160 

XINDEX=XMIN 

GOSUB  1410 

IF  L%=1  THEN  GOSUB 

IF  L%  <>  NOF%  THEN 

NEXT  L% 
ON  KEY(l)  GOSUB  1006 
ON  KEY (2)  GOSUB  1009 
KEY (2)  OFF:  KEY(l)  ON:  KEY (2) 
IF  FLAG  =  1  THEN  GOTO  1100 
GOTO  1001 


1120 

PTR%=PGSTK% (STKPTR%) 


ON 


23,1: 


PRINT  "HI" 


LOCATE 
RETURN 
FLAG  =  1 
RETURN 
FLAG  =  0 
RETURN 

i pL0T  ZER0  LINE- 

IF  (YMAX*YMIN  <0)  THEN  TEMP= (-1*YMIN/ (YMAX-YMIN) )* 150 

ELSE  GOTO  1150 

DELTAY%=CINT (TEMP) 

LINE (ORIGINX% , TEMP4%-DELTAY% ) - (ORIGINX%+7*INTERVALX% , 

TEMP4%-DELTAY%) 

IF  (XMAX*XMIN  <0)  THEN  TEMP= (-1*XMIN/ (XMAX-XMIN) ) *550 

ELSE  GOTO  1190 

DELTAX%=CINT (TEMP) 

LINE ( ORIGINX%+DELTAX% , ORIGINY% ) - ( ORIGINX%+DELTAX% , TEM 

P4%) 

i PRINT  INDECIES 

YINDEX=YMAX 
TEMP  =  (YMAX-YMIN) /5 
P%=LEN(HEAD$) :  P%=HDPOS%- (P%* . 5) 
LOCATE  1,1:  PRINT  TAB(P%)  HEAD$ 
FOR  1%  =  1  TO  6 

PRINT  USING  "###.###" ;YINDEX 

IF  (I%=  3)  THEN  FOR  J%  =  1  TO  2 : 
GOTO  12  60 

IF  (I%<>6)  THEN  FOR  J%  =  1  TO  3 : 

YINDEX=YINDEX-TEMP 

NEXT  1% 
XINDEX  =  XMIN:  XINCR= (XMAX-XMIN) /7 
LOCATE  22,1 

PRINT  USING  "########.#"; XINDEX; 
FOR  1%  =  1  TO  7 

XINDEX  =  XINDEX  +  XINCR 

PRINT  USING  "########.#"; XINDEX; 
P%=LEN(XTITLE$) :  P%=HDPOS%- (P%* . 5) 
LOCATE  23,1:  PRINT  TAB(P%)  XTITLE$ 


PRINT:  NEXT  J* 


PRINT:  NEXT  J'- 


NEXT  r 
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1350  T1%=LEN(YTITLE$) :  P%=11- (Tl%* . 5) 

1360  FOR  1%  =  1  TO  Tl% 

1370      T1$=MID$(YTITLE$,I%,1) :  LOCATE  P%+I%,1:  PRINT  Tl$ 

1380      NEXT 

1395  RETURN 

140o  ' PLOT  POINTS  ON  GRAPH 

1410  TEMP2=1/(YMAX-YMIN) :  TEMP3=1/ (XMAX-XMIN) 

142  0  IF  MRK$=""  THEN  MRK$="MARK    " 

1430  REP2$=LEFT$(MRK$,1) 

1440  T%=L%:  REP$=LEFT$ (CURVE$ (L%) , 1) 

1450  IF  T%>3  THEN  T%=T%-3  :  GOTO  1450 

1460  IF  REP$="C"  THEN  GOSUB  1870 

1470  COUNT%=0:  NDP%=0 

148  0  WHILE  (PTR%<=NOD%)  AND  (X(PTR%)  <=  XMAX) 

1490      IF  X(PTR%)  <  XMIN  THEN  GOTO  1650 

1500      OLDY%=Y%:  OLDX%=X% 

1510      NDP%  =  NDP%+1 

1520      TEMP1=ABS( (X (PTR%) -XMIN) *TEMP3) *550 

153  0      TEMP=ABS( (Y (PTR%) -YMIN) *TEMP2) *150 

1540      DELTAY%=CINT(TEMP) :  DELTAX%=CINT (TEMPI) 

1560      Y%=TEMP4%-DELTAY%  :  X%=ORIGINX%+DELTAX% 

1570      IF  COUNT%=0  THEN  OLDY%=Y% :  OLDX%=X% :  PSET(X%,Y%): 

FIRSTX%=X%:  FIRSTY%=Y%:  COUNT%=l 
1580      IF  REP$="P"  THEN  1600 
1590      LINE-(X%,Y%) ,1 
1600      IF  REP2$  <>  "M"  AND  REP  <>  "P"  THEN  GOTO  163  0 

1610  IF  T%=1  THEN  CIRCLE (X% , Y%) , 2 , 1, 0 , 6 . 28 , 5/12 

1611  IF  T%=2  THEN  LINE (X% , Y%) - (X%+6 , Y%) :  LINE- (X%+3 , Y% 
-3) :  LINE-(X%,Y%) 

1612  ELSE  LINE(X%,Y%)-(X%+6,Y%-3) ,1,B  :  PSET(X%,Y%) 
1620      GOTO  1640 

1630      PSET(X%,Y%) 

1640      IF  REP$="C"  THEN  GOSUB  2140 

1650      PTR%=PTR%+1 

1660      WEND 

172  0  RETURN 

1860  ' CUBIC  INTERPOLATION 

1870  DIM  H(200) ,A(200) ,L(200) ,U(200) ,Z(200) 

1880  FOR  1%  =  1  TO  NOD%-l 

1890      H(I%)=X(I%+1)-X(I%) 

1900      NEXT 

1910  FOR  I%=2  TO  NOD%-l 

1920      T1=Y(I%+1) *H(I%-1) 

1930      T2=Y(I%)*(X(I%+1)-X(I%-1)  ) 

1940      T3=Y(I%-1)*H(I%) 

1950      A(I%)=3*(T1-T2+T3)/(H(I%-1) *H(I%) ) 

1960      NEXT  1% 

1970  L(l)=l:  U(1)=0:  Z(1)=0 

1980  FOR  I%=2  TO  NOD% 

1990      L(I%)=2*(X(I%+1)-X(I%-1) )-H(I%-l) *U(I%-1) 

Figure  A. 15   PLOTTER. BAS  (cont.) 
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2000      U(I%)=H(I%)/L(I%) 

2010      Z(I%)=(A(I%)-H(I%-1)*Z(I%-1))/L(I%) 

2020      NEXT  1% 

2030  L(N0D%)=1:  Z(NOD%)=0:  C(NOD%)=0 

2040  FOR  J%=N0D%-1  TO  1  STEP  -1 

2050      C(J%)=Z(J%)-U(J%)*C(J%+1) 

2060      T1=(Y(J%+1)-Y(J%) )/H(J%) 

2070      T2=H(J%) *(C(J%+1)+2*C(J%) )/3 

2080      B(J%)=T1-T2 

2  090      D(J%)=(C(J%+1)-C(J%) )/(3*H(J%) ) 

2100      NEXT  J% 

2110  ERASE  H,A,L,U,Z 

2120  RETURN 

2130  ' PLOT  CUBIC 

2140  IF  PTR%+1  >  NOD%  OR  X(PTR%+1)  >  XMAX  THEN  GOTO  22  8  0 

2150  Tl=(XMAX-XMIN)/560  :  T2=X(PTR%) 

2160  TEMP1=ABS( (X(PTR%+1) -XMIN) *TEMP3 ) *550 

2170  DELTAX%=CINT (TEMPI) 

2180  XNEXT%=ORIGINX%+DELTAX% 

2190  X%=X%+2 

2  2  00  IF  X%  >=  XNEXT%  THEN  GOTO  228  0 

2210  T2=T2  +  2*T1 

2220  T4=T2-X(PTR%) 

2230  T3=Y(PTR%)+B(PTR%) *T4+C(PTR%) *T4*T4+D(PTR%) *T4*T4*T4 

224  0  TEMP=ABS( (T3-YMIN) *TEMP2) *150 

2250  DELTAY%=CINT(TEMP) :  Y%=TEMP4%-DELTAY% 

2260  IF  Y%>172  THEN  Y%=172  ELSE  IF  Y%<12  THEN  Y%=12 

2270  LINE-(X%/Y%)  :  GOTO  2190 

2280  RETURN 

2400  ' PLOT  WITH  GRIDS 

2410  FOR  1%  =  1  TO  9 

2420      Tl%=ORIGINY%+I%*15 

2430      LINE(73,Tl%)-(623/Tl%) , , ,&H4444 

2440      NEXT  1% 

2450  FOR  I%=  1  TO  13 

2460      T1%=ORIGINX%+I%*40 

2470      LINE(T1%,12)-(T1%/162) , , , &HAAAA 

2480      NEXT  1% 

2490  RETURN 

2500  END 
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NAME:   Data  Acquisition  And  Control  (DAAC) 
HEADER  for  BASICA 

FILE  NAME:   DACHDR.BAS 

DOS  DEVICE  NAME:   DAAC 

RESERVED  FUNCTION  NAMES: 

AINM,  AINS,  AINSC,  AOUM,  AOUS , 
BINM,  BINS,  BITINS,  BITOUS ,  BOUM,  BOUS , 
CINM,  CINS,  CSET,  DELAY 
RESERVED  DEF  SEG  VALUE  NAME:   DSEG 

NAMES  DEFINED  AND  USED  BY  HEADER: 

ADAPT%,  AI,  COUNT,  FOUND% , 
HNAME$,  SG%,  STAT% 


When  using  the  BASICA  Interpreter,  this  header 
must  be  executed  before  any  function  calls  are 
made  that  access  the  DAAC  adapter.   It  initializes 
a  number  of  variables  for  each  function  call.   These 
variables  are  reserved  and  should  not  be  used  except 
to  access  the  DAAC  adapter.   This  routine  also  does  a 
DEF  SEG  to  the  segment  where  the  DAAC  Device  Driver 
(DAC.COM)  is  loaded.   If  you  execute  a  DEF  SEG  to 
access  other  hardware,  you  must  DEF  SEG  to  the  segment 
of  the  DAAC  Device  Driver  before  any  subsequent 
calls  to  access  the  DAAC  adapter. 


FOUND%  =0 
SG%  =  &H2E 

'Start  searching  the  interrupt  vectors  until  you  find 
'one  that  points  to  the  DAAC  device  driver. 
'Do  a  DEF  SEG  to  that  segment. 
WHILE  ((SG%  <=  &H3E)  AND  (FOUND%  =  0)) 

DEF  SEG  =  0 

DSEG  =  PEEK(SG%)  +  PEEK(SG%  +  1)  *  256 

DEF  SEG  =  DSEG 

HNAME$="" 

FOR  AI=10  TO  17 

HNAME$  =  HNAME$  +  CHR$ (PEEK (AI) ) 

NEXT  AI 


IF  HNAME$  =  "DAAC 
0  THEN  FOUND%  =  1 
SG%  =  SG%  +  4 


AND  PEEK (18)  +  PEEK (19)  <> 


WEND 

IF  FOUND%  = 

NOT  FOUND"  : 


0  THEN  PRINT  "ERROR: 
END 


DEVICE  DRIVER  DAC.COM 


Figure  A. 16   BALCAL.BAS  -  Balance  Calibration  Program 
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'Now  initialize  all  function  name 
•to  access  the  device  driver. 


variables  for  calls 


AINM  -  PEEK(&H13)  *  256  + 

AINS  -  PEEK(&H15)  *  256  + 

AINSC  =  PEEK(&H17)  *  256  + 

AOUM  =  PEEK(&H19)  *  256  + 

AOUS  -  PEEK(SHIB)  *  256  + 

BINM  -  PEEK(&H1D)  *  256  + 

BINS  =  PEEK(SHIF)  *  256  + 

BITINS  =  PEEK(&H21)  *  256  + 

BITOUS  =  PEEK(&H23)  *  256  + 

BOUM  =  PEEK(&H25)  *  256  + 

BOUS  =  PEEK(&H27)  *  256  + 

CINM  =  PEEK(&H29)  *  256  + 

CINS  =  PEEK(&H2B)  *  256  + 

CSET  =  PEEK(&H2D)  *  256  + 

DELAY  =  PEEK(&H2F)  *  256  + 

'Finally,  execute  any  call  to 

'device  driver  from  any  former  invocation  of  BASIC. 

ADAPT %  =  0 

COUNT  =  1 

STAT%  =  0 

CALL  DELAY  (ADAPT%,  COUNT,  STAT%) 
i 

•End  of  DAAC  BASICA  Header 

i 

REM BALCAL.BAS   (CALIBRATE  BALANCE) 

REM 

DIM  L(100) ,D(100) ,Y(100) ,PM(100) ,YM(100) ,RM(100) , LBS ( 
100) 
DIM  DAT (39! 


PEEK(&H12) 
PEEK(&H14) 
PEEK(&H16) 
PEEK(&H18) 
PEEK(&H1A) 
PEEK(SHIC) 
PEEK(&H1E) 
PEEK(&H2  0) 
PEEK(&H22) 
PEEK(&H24) 
PEEK(&H2  6) 
PEEK(&H2  8) 
PEEK(&H2A) 
PEEK(&H2C) 
PEEK(&H2E) 
re-initialize  the 


=.  1 


FOR  I 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 

IF 
REM 
REM 

GOSUB  138 
NEXT  I 
CLS 
LOCATE  10 


9) ,DAT%(3  99) ,DAT1(399) ,DAT1%(3  99) 
TO  12 


I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 
I 


THEN 
THEN 
THEN 
THEN 
THEN 
THEN 
THEN 

8  THEN 

9  THEN 

10  THEN 

11  THEN 

12  THEN 


1 
2 
3 

4 
5 
6 
7 


FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$  = 

FILE$ 

FILE$ 

FILE$ 


"LIFTP" 

"DRAGP" 

"YAWP" 

"PITCHMP" 

"YAWMP" 

"ROLLMP" 

"LIFTN" 

"DRAGN" 

"YAWN" 
=  "PITCHMN" 
=  "YAWMN" 
=  "ROLLMN" 


'RECORD  CALIBRATION  DATA 


, 5: PRINT  "CALCULATING  CALIBRATION  CONSTANTS, 
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PLEASE  WAIT" 
960   GOSUB  2100   'CALCULATE  K1&K2  FOR  ALL  PRIME  GAGES 
970   GOSUB  2  640   'CONVERT  COUNTS  TO  FORCES 
980   GOSUB  3180   'CALCULATE  INTERACTION  CONSTANTS 
990   REM 

1000   REM   STORE  CALIBRATION  CONSTANTS 
1010   OPEN  "C: CONST"  FOR  OUTPUT  AS  #1 

102  0   WRITE  # 1 , INCALL , INCALD , INCALY , INCALPM , INCALYM , INCALRM 
1040   WRITE  #1,K1LP0S,K2LP0S,K1DP0S,K2DP0S,K1YP0S,K2YP0S 
1050   WRITE  #1,K1PMP0S,K2PMP0S,K1YMP0S,K2YMP0S,K1RMP0S, 

K2RMPOS 
1060   WRITE  #1,K1LNEG,K2LNEG,K1DNEG,K2DNEG,K1YNEG,K2YNEG 
1070   WRITE  #1,K1PMNEG,K2PMNEG,K1YMNEG,K2YMNEG,K1RMNEG, 

K2RMNEG 
1080   WRITE  #1,DDDL1P,DDDL2P,DYDL1P,DYDL2P,DPMDL1P,DPMDL2P 

,DYMDL1P,DYMDL2P 
109  0   WRITE  #1,DRMDL1P,DRMDL2P 
1100   WRITE  #1,DLDD1P,DLDD2P,DYDD1P,DYDD2P,DPMDD1P,DPMDD2P 

,DYMDD1P,DYMDD2P 
1110   WRITE  #1,DRMDD1P,DRMDD2P 
1120   WRITE  #1,DLDY1P,DLDY2P,DDDY1P,DDDY2P,DPMDY1P,DPMDY2P 

,DYMDY1P,DYMDY2P 
1130   WRITE  #1,DRMDY1P,DRMDY2P 
1140   WRITE  #1,DLDPM1P,DLDPM2P,DDDPM1P,DDDPM2P,DYDPM1P,DYD 

PM2  P , DYMDPM1P , DYMDPM2  P 
1150   WRITE  #1,DRMDPM1P,DRMDPM2P 
1160   WRITE  #1,DLDYM1P,DLDYM2P,DDDYM1P,DDDYM2P,DYDYM1P,DYD 

YM2 P , DPMDYM1P , DPMDYM2 P 
117  0   WRITE  #1,DRMDYM1P,DRMDYM2P 
1180   WRITE  #1,DLDRM1P,DLDRM2P,DDDRM1P,DDDRM2P,DYDRM1P,DYD 

RM2  P , DPMDRM1P , DPMDRM2  P 
119  0   WRITE  #1,DYMDRM1P,DYMDRM2P 
1200   WRITE  #1,DDDL1N,DDDL2N,DYDL1N,DYDL2N,DPMDL1N,DPMDL2N 

, DYMDL1N , DYMDL2N 
1210   WRITE  #1,DRMDL1N/DRMDL2N 
1220   WRITE  #1,DLDD1N,DLDD2N/DYDD1N,DYDD2N,DPMDD1N/DPMDD2N 

,DYMDD1N,DYMDD2N 
123  0   WRITE  #1,DRMDD1N,DRMDD2N 
1240   WRITE  #1/DLDY1N,DLDY2N,DDDY1N,DDDY2N,DPMDY1N,DPMDY2N 

,DYMDY1N,DYMDY2N 
1250   WRITE  #1,DRMDY1N,DRMDY2N 

12  60   WRITE  #1,DLDPM1N,DLDPM2N,DDDPM1N,DDDPM2N,DYDPM1N,DYD 

PM2N,DYMDPM1N 

127  0   WRITE  #1,DYMDPM2N,DRMDPM1N,DRMDPM2N 

128  0   WRITE  #1,DLDYM1N,DLDYM2N,DDDYM1N,DDDYM2N,DYDYM1N,DYD 

YM2N,DPMDYM1N 
1290   WRITE  #1,DPMDYM2N,DRMDYM1N,DRMDYM2N 

13  0  0   WRITE  #1,DLDRM1N,DLDRM2N,DDDRM1N,DDDRM2N/DYDRM1N,DYD 

RM2N,DPMDRM1N 
1310   WRITE  #1,DPMDRM2N,DYMDRM1N,DYMDRM2N 
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1320   CLOSE  #1 

13  3  0   CLS:  LOCATE  10,5 

1340   PRINT  "  CALIBRATION  COMPLETE" 

1350   END 

13  60   REM 

137  0   REM 

13  80  REM RECORD  CALIBRATION  DATA •— 

1401  COLOR  15,1:  KEY  OFF:  CLS 

1402  PRINT  "CALIBRATION  FOR  THE  LOADING  OF  THE  ";FILE$;" 
COMPONENT" 

1403  PRINT:  INPUT  "AFTER  ALL  AMPLIFIERS  ARE  ZEROED  PRESS 
RETURN" ;X 

1404  CLS 

1405  PRINT  "   ZD        ZL       ZY         ZPM       ZYM 
ZRM  " 

1406  STAT%=0:   MODE%=0 :  STOR%=0 :  COUNT=100:  RATE=500 

1407  ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=3 

14  08  CALL  AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

1409  ZD  =  0:ZPM=0:ZL=0:ZYM=0 

1410  IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 
STAT%:END 

1411  FOR  J  =  0  TO  396  STEP  4 

1412  DAT(J)=(DAT%(J)/204.8)-10 

1413  ZD  =ZD  +  DAT (J) 

1414  NEXT  J 

1415  ZD  =  ZD/100 

1416  FOR  J  =  1  TO  397  STEP  4 

1417  DAT(J)=(DAT%(J)/204.8)-10 

1418  ZL  =  ZL  +  DAT (J) 

1419  NEXT  J 

1420  FOR  J  =  2  TO  398  STEP  4 

1421  DAT(J)=(DAT%(J)/204.8)-10 

1422  ZPM  =  ZPM  +  DAT (J) 
142  3  NEXT  J 

1424  FOR  J  =  3  TO  399  STEP  4 

1425  DAT(J)=(DAT%(J)/204.8)-10 

1426  ZYM  =  ZYM  +  DAT (J) 
14  27  NEXT  J 

1428  STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

1429  ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=1 
14  3  0  CALL  AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

1431  ZY  =0:ZRM=0 

1432  IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 
STAT% : END 

1435  FOR  J  =  0  TO  198  STEP  2 

1436  DATl(J)=(DATl%(J)/204.8)-10 

1437  ZRM  =  ZRM  +  DAT1 ( J) 
14  3  8  NEXT  J 
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ZRM  =  ZRM/100 

FOR  J  =  1  TO  199  STEP  2 

DATl(J)=(DATl%(J)/204.8)-10 

ZY  =  ZY  +  DAT1(J) 

NEXT  J 

ZL=ZL/100 : ZPM=ZPM/100 : ZYM=ZYM/100 : ZY=ZY/100 


LOCATE  3,1:  PRINT  USING 
NT  USING  "+#.###" ;ZL 


••  +  #.###■■  ;ZD:    LOCATE    3,10:    PRI 


3,19 
3,28 


3,37 
3,46 


PRINT 
PRINT 
PRINT 
PRINT 


USING 
USING 
USING 
USING 


"  +  # 
"  +  # 
"  +  # 
"  +  # 


###";zy 
###";zpm 


###" 


ZYM 
ZRM 


PLACING  ALL  CAL  SWITCHES  TO  +  SETTING  HIT 


CALD 


CLL 


CALY 


CALPM 


CALYM 


MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

g.—     a,      r>u*\iT  r\<2.—     r\  •      itp-ot  9 —     q.      CHANHI%  =  3 


LOCATE 

LOCATE 

LOCATE 

LOCATE 

REM 

PRINT : PRINT 

INPUT  "AFTER 

RETURN" ;X 

PRINT 

PRINT  " 

CALRM" 

STAT%=0 

ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%= 

CALL  AINSC (ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

CALD=0 : CLL=0 : CALYM=0 : CALPM=0 

IF  STAT%<>  0,THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT% : END 

FOR  J  =  0  TO  396  STEP  4 

DAT(J)=(DAT%(J)/2  04.8)-10 

CALD  =CALD  +  DAT (J) 

NEXT  J 

CALD  =  CALD/ 100 

FOR  J  =  1  TO  397  STEP  4 

DAT(J)=(DAT%(J)/2  04.8)-10 

CLL  =  CLL  +  DAT (J) 

NEXT  J 

FOR  J  =  2  TO  3  98  STEP  4 

DAT(J)=(DAT%(J)/2  04.8)-10 

CALPM  =  CALPM  +  DAT (J) 

NEXT  J 

FOR  J  =  3  TO  3  99  STEP  4 

DAT(J)=(DAT%(J)/2  04.8)-10 

CALYM  =  CALYM  +  DAT (J) 

NEXT  J 

STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=1 

CALL  AINSC ( ADAPT% , DEVICE% , CHANLO% , CHANHI% , CTRL% , MODE% 

, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

CALY  =0:CALRM=0 

IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 

STAT%:END 
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FOR  J  =  0  TO  198  STEP  2 

DAT1 (J) = (DAT1% ( J) /204 . 8) -10 

CALRM  =  CALRM  +  DAT1(J) 

NEXT  J 

CALRM  =  CALRM/ 100 

FOR  J  =  1  TO  199  STEP  2 

DATl(J)=(DATl%(J)/2  04.8)-10 

CALY  =  CALY  +  DAT1(J) 

NEXT  J 

CLL=CLL/ 100: CALPM=CALPM/ 100: CALYM=CALYM/ 100: CALY=CALY 

/100 

PRINT  USING  "+#.###" ;CALD 
»+#.###" ;CLL 


LOCATE  10,1: 
LOCATE  10,10: 
LOCATE  10,19 
LOCATE  10,28 
LOCATE  10,37 
LOCATE  10,46 
REM 
REM 


PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  USING 


'+#.###" ;CALY 
'+#.###" ;CALPM 
'+#.###" ;CALYM 
»+#.###", -CALRM 


IF 
IF 
IF 
IF 
IF 
IF 


THEN  INCALL  =  CLL  -  ZL 

THEN  INCALD  =  CALD  -  ZD 

THEN  INCALY  =  CALT  -  ZY 

THEN  INCALPM  =  CALPM  -  ZPM 

THEN  INCALYM  =  CALYM  -  ZYM 

THEN  INCALRM  =  CALRM  -  ZRM 
PRINT: INPUT  "PRESS  ENTER  (RETURN)  TO  CONTINUE" ;X 
CLS: PRINT  "RETURN  ALL  CAL  SWITCHES  TO  CENTER  POSITION" 
PRINT : 

PRINT  "LOAD  THE  ";FILE$;"  GAGE  FROM  0  TO  2  0  POUNDS  AND 
THEN  BACK  TO  0" 

PRINT  "POUNDS  IN  1  POUND  INCREMENTS.  THERE  SHOULD  BE 
TWO  READINGS  FOR" 

PRINT  "EACH  POUND  WEIGHT  EXCEPT  FOR  THE  LOAD  AT  20 
POUNDS  WHICH  WILL  HAVE" 
PRINT  "ONLY  ONE  READING" 
PRINT 

PRINT  "PRESS  THE  Fl  KEY  WHEN  THE  LOADING  IS  FINISHED" 
PRINT  "PRESS  THE  F2  KEY  WHEN  READY  TO  RECORD  THE  DATA 
FOR  THAT  LOAD" 

PRINT:  INPUT  "PRESS  ENTER  (RETURN) 
CLS 

YAW 


TO  CONTINUE ";X 


PRINT  "  DRAG 
ROLL   " 
PRINT  " 
MOM.  " 
SOAP  =  0: 
FOR  K  =  1 
ON  KEY(l) 
ON  KEY (2) 


LIFT 


N=4 

TO  100 
GOSUB  1830 
GOSUB  1841 


PITCH 


MOM. 


YAW 


MOM. 


'SET  STOP  FLAG 
'RECORD  DATA 


KEY(l)  ON:  KEY (2)  ON 
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1760  IF  SOAP  =  2  THEN  GOTO  1780 

1770  IF  SOAP  =  1  THEN  GOTO  192  0 

1775  GOTO  1730 

1780  SOAP  =  0 

1790  NOD%  =  K 

1800  NEXT  K 

1810  GOTO  1920 

183  0  REM   SET  STOP  FLAG 

1838  SOAP  =  1 

1839  RETURN 

184  0  REM  STEPS  TO  RECORD  DATA  ' 

1841  STAT%=0:   MODE%=0:  STOR%=0:  COUNT=100:  RATE=500 

1842  ADAPT%=  0:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=3 

1843  CALL  AINSC(ADAPT%/DEVICE%,CHANLO%/CHANHI%/CTRL%,MODE% 
, STOR% , COUNT , RATE , DAT% ( 0 ) , STAT% ) 

1844  D(K)=0:L(K)=0:YM(K)=0:PM(K)=0 

1845  IF  STAT%<>  0,THEN  PRINT  USING  "EXECUTION  ERROR  ###" ; 
STAT% : END 

1846  FOR  J  =  0  TO  396  STEP  4 

1847  DAT(J)=(DAT%(J)/204.8)-10 

1848  D(K)  =D(K)  +  DAT(J) 

1849  NEXT  J 

1850  D(K)  =  D(K)/100 

1851  FOR  J  =  1  TO  397  STEP  4 

1852  DAT(J)=(DAT%(J)/204.8)-10 

1853  L(K)  =  L(K)  +  DAT(J) ' 

1854  NEXT  J 

1855  FOR  J  =  2  TO  398  STEP  4 

1856  DAT(J)=(DAT%(J)/204.8)-10 

1857  PM(K)  =  PM(K)  +  DAT (J) 

1858  NEXT  J 

1859  FOR  J  =  3  TO  399  STEP  4 

1860  DAT(J)=(DAT%(J)/204.8)-10 

1861  YM(K)  =  YM(K)  +  DAT(J) 
18  62  NEXT  J 

1863  STAT%=0:   MODE%=0:  STOR%=0 :  COUNT=100:  RATE=500 
18  64  ADAPT%=  1:  DEVICE%=  9:  CHANLO%=  0:  CTRL%=  0:  CHANHI%=1 
18  65  CALL  AINSC(ADAPT%,DEVICE%,CHANLO%,CHANHI%,CTRL%,MODE% 
, STOR% , COUNT , RATE , DAT1% ( 0 ) , STAT% ) 

1866  Y(K)  =0:RM(K)=0 

1867  IF  STAT%<>  0  THEN  PRINT  USING  "EXECUTION  ERROR  ###"; 
STAT% : END 

1868  FOR  J  =  0  TO  198  STEP  2 

1869  DATl(J)=(DATl%(J)/204.8)-10 

1870  RM(K)  =  RM(K)  +  DAT1(J) 

1871  NEXT  J 

1872  RM(K)  =  RM(K)/100 

1873  FOR  J  =  1  TO  199  STEP  2 

1874  DATl(J)=(DATl%(J)/204.8)-10 

1875  Y(K)  =  Y(K)  +  DAT1(J) 
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NEXT  J 

L(K)=L(K)/100:PM(K)=PM(K)/100:YM(K)=YM(K)/100:Y(K)=Y( 
K)/100 
IF  N>  23  THEN  N=4 : CLS : PRINT  "  DRAG      LIFT       YAW 
PITCH      YAW      ROLL  " : PRINT  " 
MOM.       MOM.     MOM.  " 
LOCATE  N,l:  PRINT  USING  "+#.###" ;D (K) 
LOCATE  N,10:  PRINT  USING  "+#.###» ;L(K) 

"+#.###";Y(K) 
"+#.###»;PM(K) 
"+#.###";YM(K) 
»+#.###»;RM(K) 


N,10 
N,19 
N,28 
N,37 
N,46 


PRINT 
PRINT 
PRINT 
PRINT 


USING 
USING 
USING 
USING 
USING 


DIFFD=(D(NOD%)-D(l) ) :DIFFY=(Y( 


LOCATE 
LOCATE 
LOCATE 
LOCATE 
N=N+1 
SOAP  =  2 

RETURN    "GO  BACK  FOR  NEXT  DATA  POINT 
'—CORRECT  FOR  DRIFT- 
DIFFL  =(L(NOD%)-L(l) ) 
NOD%)-Y(l) ) 

DIFFPM=(PM(NOD%)-PM(l) ) : DIFFYM= (YM(NOD%) -YM( 1) ) :DIFFR 
M=(RM(NOD%)-RM(l) ) 
D=NOD%-l 

CORL=DIFFL/D : CORD=DIFFD/D : CORY=DIFFY/D : CORPM=DIFFPM/D 
: CORYM=DIFFYM/D : CORRM=DIFFRM/D 
A=l 
FOR  K  =  2  TO  NOD% 

L(K)=L(K)-(A*CORL) 

D(K)=D(K)-(A*CORD) 

Y(K)=Y(K)-(A*CORY) 

PM(K) =PM (K) - (A*CORPM) 

YM(K) =YM (K) - (A*CORYM) 

RM(K) =RM (K) - (A*CORRM) 

A  =  A+l 
NEXT  K 

i — -WRITE  DATA  TO  FILE 

FILES  =  "C:"+FILE$ 

OPEN  FILE$  FOR  OUTPUT  AS  #1 

WRITE  #1,ZL, ZD,ZY,ZPM,ZYM,ZRM,NOD% 

WRITE  # 1 , CLL , CALD , CALY , CALPM , CALYM , CALRM 

IF  1=1  OR  1=2  OR  1=3  THEN  GOTO  2  079 

IF  1=7  OR  1=8  OR  1=9  THEN  GOTO  2  067 


LBS  = 

FOR  J 

IF 

IF 

IF 


0 

=  1 
1=4 
1=5 
1=6 


TO  11 

OR  1=10  THEN  Dl=20 

OR  1=11  THEN  Dl=4.5 

OR  1=12  THEN  Dl=11.5 
MOM  =  LBS*D1  'CONVERT  TO  MOMENTS 
IF  1=10  OR  1=11  OR  1=12  THEN  MOM=MOM*(-l) 
WRITE  #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,MOM 
LBS  =  LBS+1 
NEXT  J 
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1986  LBS  =  9 

1987  FOR  J  =  12  TO  NOD% 

1988  IF  1=4  OR  1=10  THEN  Dl=20 

1989  IF  1=5  OR  1=11  THEN  Dl=4 . 5 

1990  IF  1=6  OR  1=12  THEN  Dl=11.5 

1991  MOM  =  LBS*D1  'CONVERT  TO  MOMENTS 

1992  IF  1=10  OR  1=11  OR  1=12  THEN  MOM=MOM*(-l) 

1993  WRITE  #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,MOM 

1994  LBS  =  LBS-1 

1995  NEXT  J 

1996  CLOSE  #1 

1997  RETURN     'GO  BACK  TO  SET  UP  NEXT  CALIBRATION 

2067  LBS  =  0 

2068  FOR  J  =  1  TO  21 

2069  WRITE  #1, L( J) , D (J) , Y ( J) , PM(J) , YM( J) ,RM( J) , LBS 
2  07  0     LBS  =  LBS  -  1 

2  071     NEXT  J 

2072  LBS  =  -19 

2073  FOR  J  =  22  TO  NOD% 

2074  WRITE  #1, L(J) , D(J) , Y (J) , PM(J) , YM(J) ,RM(J) , LBS 
2  075     LBS  =  LBS  +  1 

2  076     NEXT  J 

2077  CLOSE  #1 

2078  RETURN     'GO  BACK  TO  SET  UP  NEXT  CALIBRATION 

2079  LBS  =  0 

2080  FOR  J  =  1  TO  21 

2081  WRITE  #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,  LBS 
2  082     LBS  =  LBS  +  1 

2  083     NEXT  J 

2084  LBS  =  19 

2  085  FOR  J  =  22  TO  NOD% 

2086     WRITE  #1, L( J) , D ( J) , Y ( J) , PM ( J) , YM ( J) , RM ( J) , LBS 

2  087     LBS  =  LBS  -  1 

2  088     NEXT  J 

2089  CLOSE  #1 

2090  RETURN     'GO  BACK  TO  SET  UP  NEXT  CALIBRATION 

2100  REM 

2110  REM  CALCULATE  K1&K2  FOR  PRIME  GAGES 

2120  FOR  I  =  1  TO  12 

213  0     IF  I  =  1  THEN  FILE$  =  "LIFTP" 

214  0     IF  I  =  2  THEN  FILE$  =  "DRAGP" 
2150     IF  I  =  3  THEN  FILE$  =  "YAWP" 
2160     IF  I  =  4  THEN  FILE$  =  "PITCHMP" 
217  0     IF  I  =  5  THEN  FILE$  =  "YAWMP" 
2180     IF  I  =  6  THEN  FILE$  =  "ROLLMP" 
2190     IF  I  =  7  THEN  FILE$  =  "LIFTN" 
2200     IF  1=8  THEN  FILE$  =  "DRAGN" 
2210     IF  I  =  9  THEN  FILE$  =  "YAWN" 
2220     IF  I  =  10  THEN  FILE$  =  "PITCHMN" 
223  0     IF  I  =  11  THEN  FILE$  =  "YAWMN" 
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IF  I  =  12  THEN  FILE$  =  "ROLLMN" 
GOSUB  2295    'READ  FILES  AND  PERFORM  CALCULATIONS 
NEXT  I 

RETURN      'GO  BACK  TO  CONVERT  COUNTS 
REM 
REM 

A=0 : B=0 : C=0 : D=0 : E=0 : 
OPEN  "C:"+FILE$  FOR  INPUT  AS  #2 
INPUT  #2 , ZL, ZD, ZY, ZPM, ZYM, ZRM,NOD% 
INPUT  #  2 , CLL , CALD , CALY , CALPM , CALYM , CALRM 
FOR  J  =  1  TO  NOD% 

INPUT  #2,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 


A 
B 
C 
D 

E 

NEXT  J 
CLOSE  #2 
K2 
Kl 


IF  I  =  1  OR  I  = 
IF  I  =  2  OR  I  = 
IF  I  =  3  OR  I  = 
IF  I  =  4  OR  I  = 
IF  I  =  5  OR  I  = 
IF  I  =  6  OR  I  = 
A  +  (XA2) 
B  +  (XA3) 
C  +  (XA4) 
D  +  (X*LBS(J) ) 
E  +  (X*X*LBS(J) ) 


7  THEN  X  =  L(J) 

8  THEN  X  =  D(J) 

9  THEN  X  =  Y(J) 

10  THEN  X  =  PM(J) 

11  THEN  X  =  YM(J) 

12  THEN  X  =  RM(J) 


IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 
IF 


((D/A)-(E/B))/((B/A)-(C/B)) 

(D/A)  -  (K2*(B/A)) 

=  1  THEN  K1LPOS  =  Kl :  K2LPOS  = 
THEN  K1DPOS  =  Kl:  K2DPOS  = 
THEN  K1YPOS  =  Kl:  K2YPOS  = 
THEN  K1PMPOS  =  Kl:  K2PMPOS 
THEN  K1YMPOS  =  Kl:  K2YMPOS 
THEN  K1RMPOS  =  Kl 
THEN  K1LNEG  =  Kl : 
THEN  K1DNEG  =  Kl: 
THEN  K1YNEG  =  Kl: 

10  THEN  K1PMNEG 

11  THEN  K1YMNEG 

12  THEN  K1RMNEG 


2 

3 
4 
5 
6 
7 
8 
9 


K2RMPOS 

K2LNEG  = 

K2DNEG  = 

K2YNEG  = 

=  Kl:  K2PMNEG  = 

=  Kl:  K2YMNEG  = 

=  Kl:  K2RMNEG  = 


K2 
K2 
K2 

=  K2 
=  K2 
=  K2 
K2 
K2 
K2 


K2 
K2 
K2 


RETURN 
REM 


'GO  BACK  TO  CALCULATE  K1&K2  FOR  NEXT  FILE 


REM  CONVERT  COUNTS  TO  FORCES 

FOR  I  =  1  TO  12 

IF  I  =  1  THEN  FILE$ 
IF  I  =  2  THEN  FILE$ 
IF  I  =  3  THEN  FILE$ 
IF  I  =  4  THEN  FILE$ 
IF  I  =  5  THEN  FILE$ 
IF  I  =  6  THEN  FILE$ 
IF  I  =  7  THEN  FILE$ 


"LIFTP" 

"DRAGP" 

"YAWP" 

"PITCHMP" 

"YAWMP" 

"ROLLMP" 

"LIFTN" 
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2740     IF  1=8  THEN  FILE$  =  "DRAGN" 

2750     IF  I  =  9  THEN  FILE$  =  "YAWN" 

2760     IF  I  =  10  THEN  FILE$  =  "PITCHMN" 

2770     IF  I  =  11  THEN  FILE$  =  "YAWMN" 

2780     IF  I  =  12  THEN  FILE$  =  "ROLLMN" 

2790  GOSUB  2835    'READ  IN  FILE  AND  CONVERT 

2800  NEXT  I 

2810  RETURN     'GO  BACK  TO  CALCULATE  INTERACTON  CONSTANTS 

2820  REM 

2830  REM 

2835  A=0:B=0:C=0:D=0:E=0:F=0 

2840  OPEN  "C:"+FILE$  FOR  INPUT  AS  #2 

2850  INPUT  #2,ZL,ZD,ZY,ZPM,ZYM,ZRM,NOD% 

2  8  60  INPUT  #2,CLL,CALD,CALY,CALPM,CALYM,CALRM 
2870  FOR  J  =  1  TO  NOD% 

2880     INPUT  #2,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 

2890     A=( (INCALL/CLL-ZL)*(L(J)-ZL) ) 

2900     B=( (INCALD/CALD-ZD) *(D(J)-ZD) ) 

2910     C=( (INCALY/CALY-ZY)*(Y(J)-ZY) ) 

2920     D=( (INCALPM/CALPM-ZPM) *(PM(J)-ZPM) ) 

293  0     E=( (INCALYM/CALYM-ZYM)*(YM(J)-ZYM) ) 

2940     F=( (INCALRM/CALRM-ZRM)*(RM(J)-ZRM) ) 

2950     IF  L(J)  <  0  THEN  Kl  =  K1LNEG:  K2  =  K2LNEG  ELSE  Kl 

=  K1LP0S:K2  =  K2LPOS 
2960     L(J)=(K1*A)  +  (K2*(AA2)) 
2970     IF  D(J)  <  0  THEN  Kl  =  K1DNEG:  K2  =  K2DNEG  ELSE  Kl 

=  K1DP0S:K2  =  K2DPOS 
2980     D(J)=(K1*B)  +  (K2*(BA2)) 
2990     IF  Y(J)  <  0  THEN  Kl  =  K1YNEG:  K2  =  K2YNEG  ELSE  Kl 

=  K1YP0S:K2  =  K2YPOS 
3000     Y(J)=(K1*C)  +  (K2*(CA2)) 
3010     IF  PM(J)  <  0  THEN  K1=K1PMNEG: K2=K2PMNEG  ELSE  K1=K1 

PMPOS:K2=K2PMPOS 
3020     PM(J)=(K1*D)  +  (K2*(DA2)) 

3  03  0     IF  YM(J)  <  0  THEN  K1=K1YMNEG: K2=K2YMNEG  ELSE  K1=K1 

YMPOS : K2=K2YMPOS 
3040     YM(J)=(K1*E)  +  (K2*(EA2)) 
3  050     IF  RM(J)  <  0  THEN  K1=K1RMNEG: K2=K2RMNEG  ELSE  K1=K1 

RMPOS : K2=K2RMPOS 
3060     RM(J)=(K1*F)  +  (K2*(FA2)) 
3  07  0     NEXT  J 
3080  CLOSE  #2 

3  090  REM   SAVE  CONVERTED  COUNTS 
3100  OPEN  "B:"+FILE$  FOR  OUTPUT  AS  #1 
3110  WRITE  #l,ZL,ZD,ZY,ZPM,ZYM,ZRM,NOD% 

312  0  WRITE  #1,CLL,CALD,CALY,CALPM,CALYM/CALRM 

313  0  FOR  J  =  1  TO  NOD% 

3140     WRITE  #1,L(J) ,D(J) ,Y(J) ,PM(J) ,YM(J) ,RM(J) ,LBS(J) 
3150     NEXT  J 
3155  CLOSE  #1 
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3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3365 


RETURN 

REM' 

REM 


'GO  BACK  TO  CALCULATE  INTERACTION  CONSTANTS 

CALCULATE  INTERACTION  CONSTANTS 
FOR  I  =  1  TO  12 

IF  I  =  1  THEN  FILE$  =  "LIFTP" 

IF  I  =  2  THEN  FILE$  =  "DRAGP" 

IF  I  =  3  THEN  FILE$  =  "YAWP" 

IF  I  =  4  THEN  FILE$  =  "PITCHMP" 

IF  I  =  5  THEN  FILE$  =  "YAWMP" 

IF  I  =  6  THEN  FILE$  =  "ROLLMP" 

IF  I  =  7  THEN  FILE$  =  "LIFTN" 

IF  I  =  8  THEN  FILE$  =  "DRAGN" 

IF  I  =  9  THEN  FILE$  =  "YAWN" 

IF  I  =  10  THEN  FILE$  =  "PITCHMN" 

IF  I  =  11  THEN  FILE$  =  "YAWMN" 

IF  I  =  12  THEN  FILE$  -  "ROLLMN" 
GOSUB  3  3  65   'READ  IN  FILES  AND  CALCULATE  CONSTANTS 
NEXT  I 

RETURN      'GO  BACK  TO  SAVE  CALIBRATION 
REM 

=0:E1=0:E2=0:E3=0:E 


DATA 


3470 

3480 

3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 


(X*Y1) 

( (XA2) *Y1) 

(X*Y2) 

( (XA2) *Y2) 

(X*Y3) 

( (XA2) *Y3) 

(X*Y4) 
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3590     E4  =  E4  +  ((XA2)*Y4) 

3600     D5  =  D5  +  (X*Y5) 

3610     E5  =  E5  +  ((XA2)*Y5) 

3  620     NEXT  J 

3630  CLOSE  #2 

3640  K12  =  ( (D1/A)-(E1/B) )/( (B/A)-(C/B) ) 

3650  Kll  =  (D1/A)-(K12*(B/A) ) 

3660  K22  =  ( (D2/A)-(E2/B) )/( (B/A)-(C/B) ) 

3670  K21  =  (D2/A)-(K22*(B/A) ) 

3680  K32  =  ((D3/A)-(E3/B) )/( (B/A)-(C/B) ) 

3690  K31  =  (D3/A)-(K32*(B/A) ) 

3700  K42  =  ((D4/A)-(E4/B) )/( (B/A)-(C/B) ) 

3710  K41  =  (D4/A)-(K42*(B/A) ) 

3720  K52  =  ( (D5/A)-(E5/B) )/( (B/A)-(C/B) ) 

3730  K51  =  (D5/A)-(K52*(B/A) ) 

3740  IF  1=1  THEN  DDDL1P=K11: DDDL2P=K12 : DYDL1P=K21: DYDL2P=K 

22:DPMDL1P=K31 
3750  IF  1=1  THEN  DPMDL2P=K32 : DYMDL1P=K41: DYMDL2P=K42 : DRMDL 

1P=K51 : DRMDL2P=K52 
3760  IF  1=2  THEN  DLDD1P=K11 : DLDD2P=K12 : DYDD1P=K21: DYDD2P=K 

22:DPMDD1P=K31 
3  77  0  IF  1=2  THEN  DPMDD2P=K3  2 : DYMDD1P=K41 : DYMDD2P=K42 : DRMDD 

1P=K51 : DRMDD2P=K52 
3780  IF  1=3  THEN  DLDY1P=K11: DLDY2P=K12 : DDDY1P=K21 : DDDY2P=K 

22:DPMDY1P=K31 
3790  IF  1=3  THEN  DPMDY2P=K32 : DYMDY1P=K41: DYMDY2P=K42 : DRMDY 

1P=K5 1 : DRMDY2  P=K5  2 
3800  IF  1=4  THEN  DLDPM1P=K11 : DLDPM2P=K12 : DDDPM1P=K21 : DDDPM 

2  P=K2  2 : DYDPM1P=K3 1 
3810  IF  1=4  THEN  DYDPM2P=K32 : DYMDPM1P=K41: DYMDPM2P=K42 : DRM 

DPM1P=K51 : DRMDPM2P=K52 
3  82  0  IF  1=5  THEN  DLDYM1P=K11 : DLDYM2P=K12 : DDDYM1P=K21 : DDDYM 

2P=K22:DYDYM1P=K31 
3  83  0  IF  1=5  THEN  DYDYM2P=K32 : DPMDYM1P=K41 : DPMDYM2P=K42 : DRM 

DYM1P=K51 : DRMDYM2P=K52 
3  84  0  IF  1=6  THEN  DLDRM1P=K11 : DLDRM2P=K12 : DDDRM1P=K21: DDDRM 

2P=K22:DYDRM1P=K31 
3  850  IF  1=6  THEN  DYDRM2P=K32 : DPMDRM1P=K41 : DPMDRM2P=K4  2 : DYM 

DRM1P=K51 : DYMDRM2P=K52 
38  60  IF  1=7  THEN  DDDL1N=K11 : DDDL2N=K12 : DYDL1N=K21 : DYDL2N=K 

22:DPMDL1N=K31 
3  87  0  IF  1=7  THEN  DPMDL2N=K3  2 : DYMDL1N=K41: DYMDL2N=K4  2 : DRMDL 

1N=K51 : DRMDL2N=K52 
3880  IF  1=8  THEN  DLDD1N=K11 : DLDD2N=K12 : DYDD1N=K21 : DYDD2N=K 

22:DPMDD1N=K31 
3890  IF  1=8  THEN  DPMDD2N=K3  2 : DYMDD1N=K41 : DYMDD2N=K4  2 : DRMDD 

1N=K51 : DRMDD2N=K52 
3900  IF  1=9  THEN  DLDY1N=K11 : DLDY2N=K12 : DDDY1N=K21 : DDDY2N=K 

22:DPMDY1N=K31 
3910  IF  1=9  THEN  DPMDY2N=K3  2 : DYMDY1N=K41 : DYMDY2N=K4  2 : DRMDY 
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1N=K51 : DRMDY2N=K52 
3  92  0  IF  1=10  THEN  DLDPM1N=K11: DLDPM2N=K12 : DDDPM1N=K21: DDDP 

M2N=K22 : DYDPM1N=K31 
3  93  0  IF  1=10  THEN  DYDPM2N=K3  2 : DYMDPM1N=K41 : DYMDPM2N=K42 : DR 

MDPM1N=K51 
3  940  IF  1=10  THEN  DRMDPM2N=K52 

3  950  IF  1=11  THEN  DLDYM1N=K11: DLDYM2N=K12 : DDDYM1N=K21 : DDDY 

M2N=K2  2 : DYDYM1N=K3 1 
3960  IF  1=11  THEN  DYDYM2N=K32 : DPMDYM1N=K41 : DPMDYM2N=K4  2 : DR 

MDYM1N=K51 
3970  IF  1=11  THEN  DRMDYM2N=K52 
3980  IF  1=12  THEN  DLDRM1N=K11: DLDRM2N=K12 : DDDRM1N=K21 : DDDR 

M2N=K2  2 : DYDRM1N=K3 1 
3990  IF  1=12  THEN  DYDRM2N=K32 : DPMDRM1N=K41 : DPMDRM2N=K4  2 : DY 

MDRM1N=K51 
4000  IF  1=12  THEN  DYMDRM2N=K52 

4  010  RETURN    'GOBACK  TO  SAVE  CALIBRATION  DATA 
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APPENDIX  B 
TABLES 


INITIAL  SET-UP  FOR  THE  MODEL  8255  TRANSDUCER  AMPLIFIER 


Amp.  #    Component   Gain  Set.   Filter  Set.    Ex.  Volt 


1 

DRAG 

MAX  VAR. 

1 

2 

LIFT 

MAX  VAR. 

1 

3 

PITCH  M. 

IK 

1 

4 

YAW  M. 

MAX  VAR. 

1 

5 

ROLL  M. 

MAX  VAR. 

1 

6 

YAW 

MAX  VAR. 

1 

7 

AOA 

IK 

1 

+  5.0 
+  6.5 
+  5.0 
+  5.0 
+  5.0 
+  5.0 
+  5.0 


Table  B.l   Amplifier  Set-Up 
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Table  B.13   Loading  for  Negative  Rolling  Moment 
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